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Preface 


In  order  to  make  our  two  year  association  with  the  Air  Force  a 
memorable  one,  it  was  decided  that  any  thesis  we  did  would  involve  the 
zero  gravity  aircraft  located  at  Wright-Patterson  AFB.  With  this  as 
a  personal  goal,  we  set  out  in  search  of  a  worthy  topic.  Based  upon 
the  guidance  of  Captain  S.  W.  Johnson,  the  topic  contained  in  this 
report  was  chosen.  For  his  patience  and  assistance  in  our  time  of 
decision  we  are  most  grateful. 

The  idea  of  detonating  an  explosive  on  board  an  aircraft  in 
flight  was  admittedly  not  an  easy  one  to  propose  to  the  Air  Force. 
Through  the  offorts  of  Mr.  Donald  Griggs  of  the  Zero  Gravity  Section, 
the  obstacles  involving  the  flight  portion  of  the  investigation  were 
overcome.  Without  his  persistence  and  ability  to  overcome  the  seem¬ 
ingly  impossible,  this  thesis  would  literally  not  have  gotten  off  the 
ground. 

The  report  attempts  to  present  a  complex  phenomenon  in  a  funda¬ 
mental  manner.  This  was  done  in  order  to  make  the  principles  de¬ 
veloped  applicable  to  a  wide  range  of  explosives  and  mediums  even 
though  the  experiment  dealt  with  small  scale  explosives  and  a  single 
medium.  Since  this  wa%  to  our  knowledge,  the  first  experiment  ever 
attempted  on  cratering  in  varied  gravity  fields,  we  feel  that  it 
marks  an  important  step  into  an  hitherto  unexplored  area. 

The  assistance  that  various  individuals  gave  us  throughout  the 
diiferent.  phases  of  the  investigation  was  outstanding.  We  are  indebted 
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to  Messrs.  Jim  Taylor,  Jack  Warwick,  and  Robert  Price  for  their  en¬ 
thusiastic  support  in  helping  us  solve  our  photography  problems.  The 
assistance  given  to  us  by  Lieutenant  Larry  Pierce  of  the  Zero  Gravity 
Section  was  an  important  factor  in  +he  successful  completion  of  the 
flight  phase  of  our  investigation.  We  would  also  like  to  thank 
Captain  S.  Alexander  and  Dr.  D.  0.  Norris  for  their  assistance, 
interest,  and  support  during  the  analysis  phase  of  the  investigation. 
Finally,  we  wish  to  thank  our  wives  for  the  encouragement  and  assist¬ 
ance  we  received  throughout  the  investigation. 

Leon  K.  Moraski 

David  J.  Teal 
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Abstract 

An  experimental  investigation  was  conducted  on  the  effects  of 
gravity  on  explosion  crater  formation.  Accompanying  the  presentation 
of  the  experimental  results  is  a  discussion  of  gravity  effects  on  the 
cratering  processes  of  compaction,  scouring,  spalling,  and  gas  accel¬ 
eration.  The  discussion  suggests  that  gravity  influences  the  shear 
resistance  of  the  medium  and  the  range  of  the  particle  throwout  and 
hence  should  influence  the  crater  formation. 

To  test  this  hypothesis,  an  experiment-  was  conducted  in  a  medium 
of  dry  silica  sand  using  a  six-grain  electrically  initiated  squib  as 
the  explosive  source.  Craters  were  formed  under  0.17  g,  O.38  g, 

1.0  g,  and  2.5  g  conditions  and  at  depths  of  burst  ranging  from  sur¬ 
face  to  five  inches.  The  statistical  t-test  was  used  to  determine  if 
the  crater  depth  and  diameter  were  affected  by  the  gravity  field. 
Crater  dimensions  versus  depth  of  burst  (DOB)  were  compared  for  the 
different  gravity  fields  to  study  the  effect  on  the  optimum  and  maxi¬ 
mum  DOB.  An  equation  expressing  crater  dimensions  as  a  decreasing 
exponential  function  of  increasing  gravity  was  assumed.  Parameters 
of  the  exponential  function  were  adjusted  to  achieve  the  best  fit  in 
the  least-squares  sense. 

The  conclusions  arrived  at  during  the  investigation  are  that  (l) 
crater  depth,  diameter,  optimum  DOB  for  depth,  and  maximum  BIB  vary 
inversely  with  gravity,  (2)  the  optimum  DOB  for  diameter  is  not  a 
function  of  gravity,  and  (3)  the  sensitivity  of  changes  in  crater 
dimensions  to  changes  in  gravity  is  a  function  of  DOB. 

xi 
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AN  INVESTIGATION  OF  THE  EFFECTS  OF  GRAVITY  ON 
CRATER  FORMATION  IN  A  COHESIONLESS  MEDIUM 

I.  Introduction 

Statement  of  Purpose  and  Scope 

The  objective  of  this  investigation  was  to  determine  the  effect  of 
gravity  on  the  formation  of  explosion  craters  in  a  cohesionless  medium. 
A  discussion  based  on  principles  of  soil  and  explosion  mechanics  is 
given.  An  experiment  was  performed  in  a  cohesionless  medium  to  test 
the  validity  of  the  hypotheses  formulated  in  the  discussion.  In  the 
experiment,  explosion  craters  were  formed  by  the  detonation  of  identi¬ 
cal  sources  of  chemical  energy  at  various  depths  of  burst  and  at  four 
different  levels  of  gravity. 

Background 

The  subject  of  cratering  as  a  means  of  large  scale  soil  and  rock 
excavation  has  been  of  interest  to  engineers  for  many  years.  Numerous 
studies  have  been  made  in  order  to  obtain  a  better  understanding  of  the 
complex  phenomenon  of  crater  formation.  These  studies  have  revealed 
that  terrestrial  craters  produced  by  subsurface  explosions  have  dimen¬ 
sions  which  depend  upon  numerous  factors,  the  most  important  being 
(1)  the  energy  of  the  explosion,  (2)  the  depth  of  the  explosion,  and 
(3)  the  physical  properties  of  the  medium  in  which  the  explosion  takes 
place. 

Permanent  manned  bases  may  be  established  on  the  Moon  in  the  near 
future.  If  these  bases  include  subsurface  installations,  there  will 
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be  a  need  for  a  reliable  and.  efficient  means  of  excavation  for  the 
placement  of  the  installation  (Ref  36:113) •  Once  a  primitive  base  is 
established,  dependence  upon  local  resources  might  require  that  some 
type  of  crude  mining  and  quarrying  operations  be  initiated.  Explo¬ 
sives,  either  chemical  or  nuclear,  can  be  used  for  these  operations. 

If  explosives  are  used  in  an  extraterrestrial  environment,  the 
effects  of  environmental  factors  such  as  gravity  and  atmosphere  on  the 
performance  of  the  explosives  must  be  known  in  order  that  the  effi¬ 
ciency  of  the  explosive  as  a  means  of  excavation  can  be  properly  eval¬ 
uated.  It  is  the  purpose  of  this  report  to  present  the  results  of  an 
investigation  of  the  effects  of  gravity  on  crater  formation. 

Current  knowledge  on  crater  formation  in  other  than  a  terrestrial 
(l  g)*  field  is  very  limited.  A  few  authorities  have  discussed  in  a 
qualitative  manner  portions  of  the  problem  and  have  arrived  at  differ¬ 
ing  conclusions  concerning  the  probable  effect  of  gravity  (Ref  39:325$ 
21:12) . 

Two  quantitative  studies  have  been  made  where  the  lunar  environ¬ 
ment  of  0.17  g  and  the  absence  of  an  atmosphere  were  considered.  These 
studies  also  resulted  in  conflicting  conclusions  (Ref  27:77;  25:31*12). 
No  experimental  work  appears  to  have  been  done  on  the  problem  of 
gravity  effects  on  crater  formation. 

Discussion  of  Explosion  Crater  Formation 

After  an  investigation  into  current  theories  of  crater  formation 

*  1  g  is  the  gravity  field  on  earth  and  is  equal  to  32.2  feet  per 
second  per  second.  Throughout  this  paper,  the  varying  gravity  fields 
are  expressed  in  terms  of  the  terrestrial  gravity  field. 
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and  size  prediction,  it  was  decided  that  a  discussion  of  the  cratering 
processes  based  on  soil  and  explosion  mechanics  principles  would  be  the 
most  feasible.  One  of  the  reasons  for  this  decision  was  that  the  pro¬ 
cesses  considered  are  operative  during  the  formation  of  craters  of  any 
size  explosion  in  any  medium.  The  other  reason  was  that  the  soil 
mechanics  approach  considers  the  shear  strength  of  the  medium  which 
the  commonly  used  hydrodynamic  model  does  not  consider.  As  is  shown 
in  this  paper,  shear  strength  is  an  important  property  in  the  consider¬ 
ation  of  gravity  effects.  After  the  discussion  of  the  processes  of 
crater  formation  under  terrestrial  conditions,  certain  hypotheses  are 
offered  as  to  the  likely  effect  of  gravity  on  each  of  these  processes 
and  the  resultant  effect  on  the  final  formation  of  the  crater. 

Experimental  Procedure 

In  order  to  test  the  hypotheses  offered  as  to  the  effect  of  grav¬ 
ity  on  crater  formation,  an  experimental  phase  was  conducted.  Small 
explosive  charges  were  detonated  in  Ottawa  sand  and  crater  depths  and 
diameters  were  measured.  Tests  were  conducted  at  1.0  g  in  the  1  bora- 
tory  and  at  0.17,  O.38,  and  2.5  g  on  a  C-13IB  aircraft.  A  statistical 
comparison  was  made  between  crater  data  taken  from  the  same  DDB^  out 
at  different  gravity  conditions  to  determine  whether  varying  the  grav¬ 
ity  field  caused  a  significant  change  in  the  crater  dimensions.  In 
addition,  a  relationship  between  crater  dimensions  and  gravity  was  as¬ 
sumed  and  the  parameters  in  the  relationship  were  adjusted  to  achieve 
the  best  fit  to  the  data  in  the  least-squares  sense. 

#  DOB  refers  to  the  depth  of  burst  of  the  explosive  beneath  the 
medium-air  interface. 
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Organization 

This  paper  is  organized  to  place  emphasis  on  the  discussion  of 
crater  formation  and  the  results  from  the  experiment.  The  detailed 
experimental  procedure,  raw  data,  methods  of  statistical  analyses,  and 
drawings  of  crater  profiles  can  be  found  in  the  appendices. 
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II .  A  Discussion  of  Gravity  Effects  on  Cratering 

Tlie  purpose  of  this  section  is  to  discuss  the  effects  of  gravity 
on  crater  formation  in  a  cohesionless  medium.  Throughout  the  dis¬ 
cussion,  a  homogeneous  medium  and  identical  explosive  sources  are  as¬ 
sumed.  A  typical  crater  is  described  and  the  terminology  used  in 
cratering  is  defined.  The  formation  of  a  crater  is  then  traced  from 
the  detonation  of  the  explosive  to  the  final  depression.  The  sources 
of  mechanical  energy  released  by  the  detonation  and  the  ensuing  cra¬ 
tering  processes  are  analyzed  to  determine  the  effect  of  gravity  on 
the  crater  depth  and  diameter.  Finally,  the  relative  importance  of 
the  cratering  processes  is  considered. 

Description  of  a  Crater 

In  order  to  discuss  the  cratering  process,  it  is  first  necessary 
to  define  the  terminology  associated  with  a  crater.  Figure  1  is  a 
cross-3ectional  profile  of  a  typical  crater.  The  material  on  the 
surface  of  the  crater  is  railed  the  fallback.  This  is  the  portion  of 
the  medium  that  was  initially  ejected  from  the  crater  but  because  of 
the  high  angle  and/or  low  velocity  of  ejection,  the  medium  particles 
were  not  thrown  clear  of  the  final  formed  crater.  The  original 
ground  material  in  contact  with  the  fallback  is  called  the  rupture 
zone.  The  medium  in  the  rupture  zone  is  completely  fragmented  and 
broken  due  to  shear  and  compressive  failure.  The  next  zone  of  in¬ 
terest  is  the  plastic  zone.  In  this  zone,  the  medium  is  permanently 
displaced  but  fractures  are  r.ot  commonly  found  here.  The  region 
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beyond  the  plastic  zone  is  called  the  elastic  zone.  In  this  zone, 
the  medium  has  not  been  permanently  deformed  by  the  cratering  process 
(Ref  h2:3hh7;  3U:l5).  The  crater  depth  and  diameter  are  defined  in 
this  paper  as  shown  in  Fig.  1. 

Sources  of  Mechanical  Energy 

Craters  of  the  type  shown  in  Fig.  1  are  formed  by  the  detonation 
of  an  explosive  or  the  hypervelocity  impact  of  a  projectile.  In  this 
discussion,  only  craters  formed  by  the  detonation  of  a  chemical  ex¬ 
plosive  source  are  considered. 

The  ignition  of  the  explosi.e  is  accompanied  by  a  great  release 
of  energy  which  is  manifested  by  high  pressures  on  the  order  of 
hundreds  of  thousands  of  atmospheres .  The  detonation  process  provides 
two  sources  of  mechanical  energy,  the  first  of  which  is  a  compressive 
shock  wave  which  causes  a  sharp  discontinuity  in  the  physical  state 
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of  the  medium.  This  shocu  wave  initially  travels  faster  than  the 
highest  elastic  wave  velocity.  As  its  energy  density  decreases,  the 
shock  wave  slows  down  and  approaches  the  compressions!  elastic  wave 
velocity.  The  second  source  of  mechanical  energy  is  provided  by  the 
pressure  produced  by  the  adiabatically  expanding  gases  near  the  source 
of  the  explosion.  These  gases  are  composed  of  (l)  the  air  in  the  voids 
of  the  medium,  (2)  the  gaseous  products  of  the  detonation,  and  (3) 
the  gases  produced  by  vaporization  of  any  moisture  in  the  medium. 

These  gases  expand  at  a  lower  velocity  than  the  shock  wave. 

Cratering  Processes 

The  sources  of  mechanical  energy  contribute  to  four  cratering 
processes  that  result  in  the  formation  of  a  crater.  These  processes 
are  (i)  compaction,  (2)  scouring,  (3)  spalling,  and  (U)  gas  accelera¬ 
tion.  Figure  2  is  a  cross-section  of  the  medium  surrounding  :he 
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explosion.  It  schematically  illustrates  the  regions  in  which  these 
processes  act  during  the  formation  of  the  crater.  Gas  acceleration 
md  spalling  take  place  primarily  in  a  cone-shaped  volume  above  the 
charge  while  compaction  and  scouring  occur  in  the  remaining  spherical 
volume  surrounding  the  charge. 

Compaction.  The  action  of  the  shock  wave  as  it  passes  through 
the  medium  is  the  primary  cause  of  compaction.  The  high  compressive 
stress  characteristic  of  the  shock  wave  and  the  pressure  discontinuity 
associated  with  the  shock  front  cause  the  medium  to  be  violently  de¬ 
formed  as  the  shock  propagates  in  a  spherically  diverging  shell  from 
the  source  of  the  explosion  (Ref  iu  :3US0) .  Because  of  the  high  initial 
stress,  the  medium  near  the  point  of  detonation  behaves  hydro dynam¬ 
ically.  This  high  stress  level  is  several  orders  cf  magnitude  higher 
than  any  restoring  forces  within  the  medium  (Ref  10:1).  As  the  initial 
stress  level  in  the  medium  caused  by  the  explosive  is  not  a  function 
of  gravity,  the  hydrodynamic  behavior  in  a  cohesionless  medium  is  in¬ 
dependent  cf  gravity. 

Due  to  energy  dissipation  resulting  from  the  spherical  expansion 
of  the  shock  wave  and  the  work  being  done  on  the  medium,  the  stress 
level  of  the  shock  wave  falls  below  the  value  necessary  for  hydro- 
dynamic  behavior  of  the  medium.  At  this  level  of  stress,  the  com¬ 
pressive  strength  of  the  .Tedium  is  exceeded  and  the  medium  is  crushed. 
This  is  the  start  of  whet  has  been  described  as  the  rupture  zone. 
Because  the  compressive  strength  of  the  medium  is  a  function  primarily 
of  molecular  bonding,  it  is  not  affected  by  gravity.  Therefore,  in 
the  first  portion  of  the  non-hydrodyr*mic  region  where  the  medium 
falls  because  its  compressive  strength  is  exceeded,  a  varying  gravity 
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field  has  no  effect . 

As  the  stress  level  further  decreases,  the  medium  in  the  rupture 
zone  undergoes  shear  failure.  Tne  shear  resistance  of  a  dry  cohesion¬ 
less  medium  can  be  expressed  by  the  equation 

s  =  p  tan  (l) 


where  s  =  shear  resistance 

p  =  effective  stress  normal  to  the  plane  of  failure 
=  angle  of  internal  friction 

The  shear  resistance  of  a  cohesionless  medium  depends  on  the  weight 
of  the  overburden.  Since  weight  is  gravity  dependent,  Eq.  (l)  can 
be  expressed  as 

S=pgHcos  i  tan  <p  (2) 

where  Q  -  ir.ass  density  of  the  medium 
g  «-  acceleration  due  to  gravity 
H  *  depth  of  overburden 

i  »  angle  shear  plane  makes  with  the  horizontal 
This  last  bquation  shows  that  (l)  the  shear  resistance  is  directly 
proportional  to  the  gravity  field  and  (2)  the  amount  of  overburden 
required  to  develop  a  certain  value  of  shear  resistance  is  inversely 
proportional  to  the  gravity  field. 

In  the  rupture  zone,  sh^ur  failure  results  from  (1)  the  dynamic 
compaction  of  the  medium  in  the  region  beneath  the  detonation  and  (2) 
the  pressure  on  the  medium  in  a  lateral  direction  from  the  detonation. 
In  considering  the  dynamic  compaction  of  the  medium,  an  analogy  of 
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the  mode  of  failure  can  be  made  to  static  compaction  because  the 
strength  of  the  medium  under  consideration  is  not  appreciably  altered 
by  differing  load  rates  (Ref  61:6?).  When  a  compressive  load  acts  on 
a  dry  cohesionless  medium,  an  inverted  triangular  prism  of  the  medium 
is  forced  downward.  As  the  prism  moves,  the  medium  at  the  sides  is 
forced  laterally  outward  as  shear  failure  occurs  along  surfaces  of 
sliding  (Ref  50:355)*  This  type  of  failure  is  shown  in  Fig.  3*  Al¬ 
though  the  spherical  loading  of  the  medium  by  the  shock  front  results 


established  that  similar  medium  failure  does  occur  and  is  instrumental, 
in  the  formation  of  a  crater.  In  a  series  of  small  scale  tests  con¬ 
ducted  at  shallow  depths  of  burst,  Baldwin  noticed  that  the  crater  rims 
were  produced  mainly  by  compressional  effects  which  squeezed  the  lower 
portions  of  the  medium  up  to  the  surface  (Ref  5:121). 


Mech/GSF-65-3b 


In  bearing  capacity  tests,  a  dense  cohesionless  soil  fails  in  the 
manner  depicted  in.  Fig.  3  (Ref  53 : 168) .  By  means  of  a  series  of  these 
tests  conducted  in  various  gravity  fields,  Halajian  showed  that  the 
bearing  capacity  of  a  dry  cohesionless  soil  varies  directly  with  grav¬ 
ity.  This  result  was  attributed  to  the  shear  resistance  of  the  soil 
which  varies  inversely  with  gravity.  Based  upon  Eq.  (2)  and  the  above 
experimental  result,  it  is  suggested  that  the  depth  and  extent  of  the 
surface  of  sliding  shown  in  Fig.  j  varies  inversely  with  gravity. 

This,  in  turn,  implies  that  the  crater  depth  and  diameter  should  vary 
inversely  w:..,h  the  gravity  field. 

The  effect  of  the  lateral  pressure  of  the  shock  wave  is  schemat¬ 
ically  represented  in  Fig.  4-  As  the  lateral  pressure  acts  on  the 
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medium,  a  portion  of  the  medium  is  displaced  along  a  surface  of  sliding 
on  which  the  shear  resistance  developed  in  the  medium  resists  the  move¬ 
ment.  The  shape  of  the  surface  of  sliding  is  determined  by  the  shear 
resistance  developed  in  the  medium.  Equation  (2)  has  shown  the  shear 
resistance  to  be  directly  proportional  to  gravity  while  the  amount  of 
overburden  required  to  develop  a  certain  value  of  shear  resistance  to 
be  invSrseiy  proportional  to  gravity.  Therefore,  the  depth  and  extent 
of  the  surfaces  of  sliding  vary  inversely  with  gravity.  This  should 
result  in  the  crater  d»ipth  and  diameter  varying  inversely  with  gravity. 

In  the  remainder  of  the  zones  formed  by  the  process  of  compaction, 
the  medium  responds  plastically  or  elastically  depending  on  the  state 
of  stress  to  which  it  is  subjected.  3ecause  of  the  negligible  amount 
of  particle  displacement  in  the  plastic  and  elastic  zones,  the  con¬ 
tribution  from  these  zones  to  the  final  crater  dimensions  is  insig¬ 
nificant  compared  to  the  contribution  from  the  rupture  and  fallback 
zones.  For  this  reason,  gravity  effects  on  these  zones  are  no.,  con¬ 
sidered. 

In  certain  types  of  explosives  known  as  deflagrating  compounds, 
the  energy  of  the  explosive  is  released  by  means  of  the  rapid  burning 
of  the  explosive  material  instead  of  a  nearly  instantaneous  decompo¬ 
sition.  This  results  in  a  slower  buildup  of  pressure  in  the  detona¬ 
tion  of  a  deflagration  compound.  As  a  result,  only  a  weak  shock  wave 
is  formed.  In  this  case,  compaction  is  also  caused  by  the  stress  gen¬ 
erated  by  the  expanding  gases .  The  main  differences  in  the  compaction 
process  are  that  the  deflagration  process  has  {l)  a  lower  stress  Jevel 
acting  on  the  medium  and  (2)  a  longer  period  of  time  in  which  the  me¬ 
chanical  energy  is  active. 
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Based  on  the  above  paragraphs,  it  is  hypothesized  that  the  process 
of  compaction  is  affected  by  a  change  in  the  gravity  field.  Although 
not  all  portions  of  the  process  are  influenced  by  the  gravity  field, 
those  that  are  affected  should  result  in  the  crater  depth  and  diameter 
varying  inversely  with  the  gravity  field . 

Scouring.  The  process  of  scouring  is  caused  by  the  turbulent 
action  of  the  gases  as  they  expand  from  the  source  of  the  explosion. 
These  gases  propagate  at  a  slower  velocity  than  the  shock  wave.  By 
means  of  an  erosive  type  of  action,  the  gases  remove  the  material 
from  the  crater  cavity  and  push  it  up  and  over  the  edge  of  the  forming 
crater.  Some  of  the  material  remains  on  the  edge  of  the  crater  and 
forms  the  lower  portion  of  the  crater  rim.  Scouring  is  normally  of 
minor  importance  compared  to  the  other  cratering  processes.  However, 
in  cases  when  the  energy  is  released  by  the  burning  of  a  deflagrating 
explosive,  a  weak  shock  is  formed  and  scouring  becomes  more  important 
in  crater  formation. 

During  scouring,  the  medium  fails  in  shear  because  the  gas  pre- 
sures  exceed  the  shear  resistance  of  the  medium.  Since  the  effective¬ 
ness  of  scouring  is  related  to  the  shear  resistance,  which,  in  turn, 
is  a  function  of  the  gravity  field,  the  effectiveness  of  scouring  at 
any  given  DOB  varies  inversely  with  the  gravity  field.  This  implies 
that  the  crater  depth  and  diameter  vary  inversely  with  the  gravity 
field. 

Spalling.  The  process  of  spalling  is  caused  by  the  interaction 
of  the  shock  wave  with  the  medium-air  interface.  To  determine  the  ef¬ 
fect  of  gravity  on  this  process,  it  is  necessary  to  trace  the 
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propagation  of  the  shock  wave  through  the  medium.  The  vertical  line 
in  Fig,  5  represents  a  plane  shock  front  traveling  with  velocity  U  in 
the  x-direction  through  the  medium.  The  undisturbed  medium  in  front 
of  the  shock  is  in  a  state  described  by  density p^  ,  particle  velocity 
l’-o>  stress  ,  and  internal  energy  E0.  The  shocked  medium  is  des¬ 
cribed  byp^  ,  u^,  j  ,  and  E^.  The  equations  that  describe  the 
medium  under  the  influence  of  a  shock  wave  are  the  Rankine-Hugoniot 
jump  conditions 


^(u-u0|  --P 


U-u 


(3) 


^(u-Ut)(u,-Lt) 


(U) 


These  equations  express  the  conservation  of  mass,  momentum,  and  energy 
across  the  shock  front.  They  relate  the  state  of  the  shocked  material 
to  the  state  of  the  undisturbed  material.  None  of  the  variables  des¬ 
cribed  in  the  Rankine-Hugoniot  jump  conditions  depend  on  the  weight 
of  the  material  concerned.  The  variables  are  described  in  terms  of 
mass  and  consequently  gravity  has  no  effect  on  them.  This  means  that 
a  shock  wave  is  independent  of  gravity  when  traveling  through  a 
medium  and  transfers  the  same  amount  of  energy  to  a  given  particle  in 
a  cohesion} ess  medium  regardless  of  the  gravity  value. 
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When  the  shock  wave  reaches 


the  medium-air  interface,  the 
boundary  condition  that  the  normal 
stress  on  the  medium  is  zero  at 
all  times  must  be  satisfied.  To 
satisfy  this  condition,  a  nega¬ 
tive  stress  wave,  or  rarefaction 
wave,  must  be  reflected  back  into 
the  medium.  This  reflected  wave 
counteracts  the  stress  produced 
by  the  incident  compressive 
wave.  This  action  is  shown 
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Figure  5 

Plane  Shock  Front 

l  From  Ref  26:66) 


graphically  in  Fig.  6.  In  this 

figure,^  is  an  exponential  stress  wave  which  represents  stress  as  a 
function  of  distance  from  the  point  of  detonation .  The  vertical  por¬ 
tion  of  the  stress  wave  is  the  shock  front.  The  wave  reaches  tha 
surface  and  is  then  reflected  as  shown  in  Fig.  6,  a,  b,  and  c.  The 
stress  level  of  the  rarefaction  is  opposite  in  sign  and  equal  in  mag¬ 
nitude  to  the  stress  level  of  the  incident  compressive  wave  at  the 
interface.  This  keeps  the  normal  stress  equal  to  aero.  4s  is  indi¬ 
cated  by  the  line  Afi  in  Fig.  6,  c,  the  algebraic  sum  of  the  stresses 
results  in  a  net  tensile  stress  in  the  medium  as  the  rarefaction 
propagates  back  toward  the  point  of  detonation.  Sines  a  cohes Ionic ss 
medium  exhibits  no  tensile  strength,  spalling  occurs  when  the  tensile 
stress  exceeds  the  restoring  force  of  gravity  which  acts  on  the  co- 
hes.ionless  particles.  The  particles  are  then  separated  from  toe 
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Figure  6 


Graphical  Illustration  ol'  Spalling  Process 
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remainder  of  the  medium  as  shown  in  Fig.  6,  d.  The  velocity  of  the 
particles,  V]_,  depends  on  the  energy  transferred  to  it  by  the  tensile 
wave  and  the  magnitude  of  the  restoring  force.  This  process  can  re¬ 
peat  itself  whenever  a  new  surface  is  created  if  the  tensile  stress, 
or  force,  exceeds  the  restoring  force  of  gravity.  The  Rankine-Hugoniot 
jump  conditions  have  shown  that  the  energy  transfer  to  the  medium  is 
independent  of  gravity.  Therefore,  a  change  in  the  gravity  field  and 
the  consequent  change  in  the  restoring  force  should  result  in  (l)  a 
change  in  the  initial  velocity  of  the  particles  which  were  originally 
spalled  and  (2)  a  change  in  the  amount  of  spalled  particles.  In  each 
case,  the  change  will  vary  inversely  with  the  change  in  gravity. 

Some  of  the  spalled  particles  are  thrown  clear  of  the  final 
crater  while  others  return  to  the  depression  as  fallback.  As  seen  in 
Fig.  1,  the  amount  of  fallback  significantly  affects  the  depth  of  the 
crater  but  not  the  diameter.  The  amount  of  material  that  becomes 
fallback  depends  on  the  range  of  the  spalled  particles  which  is  ex¬ 
pressed  by 

2 

R  -  V0  sin  2©o  (6) 

S 

where  R  «  range  of  particle  in  free  flight 

V0  -  initial  particle  velocity 

0o  -  angle  of  initial  particle  heading  as  measured  from  the 
horizontal 

g  ■  acceleration  due  to  gravity 

From  Eq.  (6)  it  is  evident  that  the  range  of  a  particle  in  free  flight 
is  inversely  proportion?!  to  gravity  if  the  remaining  tents  are  as¬ 
sumed  to  be  unchanged.  However,  it  has  been  stated  that  the  initial 
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velocity  of  the  spalled  particles  should  vary  inversely  with  a  change 
in  gravity.  This  would  also  result  in  an  additional  change  in  the 
particle  range  in  the  same  manner  (increase  or  decrease)  as  a  change 
in  the  gravity  field.  As  gravity  is  reduced,  the  resulting  increase 
in  range  results  in  a  decrease  in  the  amount  of  fallback.  Less  fall¬ 
back  results  in  a  greater  crater  depth.  In  addition,  the  increase  in 
the  amount  of  spalled  particles  that  result  from  a  decrease  in  the 
gravity  field  also  tends  to  increase  the  crater  depth.  Therefore,  it 
is  expected  that  the  crater  depth  varies  inversely  with  gravity  when 
spalling  is  considered.  Furthermore,  the  crater  diameter  is  not  ap¬ 
preciably  changed  during  the  spalling  process  when  gravity  is  varied 
as  fallback  does  not  influence  the  crater  diameter  to  a  significant 
degree. 

Gas  Acceleration.  Gas  acceleration  is  caused  by  the  escape  of 
the  expanding  gases  through  the  medium  into  the  air.  These  gases 
impart  a  force  to  the  particles  of  the  medium  and  accelerate  them  in 
accordance  with  Newton’s  second  law.  The  particles  accele.-ated  by  the 
gases  are  of  two  types.  First,  certain  spalled  particles,  because  of 
their  low  initial  velocity,  go  through  a  period  of  free  fall  in  the 
vicinity  of  the  forming  crater.  As  the  gases  escape  from  the  crater 
cavity,  they  impart  a  second  positive  acceleration  to  the  particles 
(Ref  U2i3U55)*  'Hie  final  range  of  these  particles  is  affected  by 
changes  in  the  gravity  tern  in  Eq .  (6)  and  results  in  the  depth  of  the 
crater  being  inversely  proportional  to  gravity. 

The  second  types  of  particles  accelerated  by  expanding  gases  are 
those  that  have  net  been  spalled  t  have  had  a  motion  imparted  to 
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then  by  the  passage  of  the  shock  wave.  In  this  instance  the  medium 
moves  as  one  mass  and  the  particles  remain  in  contact  with  each  other 
after  the  passage  of  the  shock  wave  (Ref  l|2:3d55)»  The  vertical 
shear  resistance  of  the  medium  now  comes  into  consideration  as  this 
resistance  must  be  overcome  before  the  particles  can  be  accelerated 
by  the  expanding  gases. 

Mediums  have  shear  resistance  on  vertical  planes  because  of  the 
lateral  pressures  within  the  material.  In  a  cohesionless  medium,  the 
lateral  pressure  acts  parallel  to  the  plane  on  which  the  vertical 
stress  actb .  This  vertical  stress  is  caused  by  the  unit  weight  of  the 
soil.  Figure  7  shows  a  stress  diagram  of  vertical  and  lateral  stress¬ 
es  acting  on  a  small  element  of  the  cohesionless  medium. 


P 

Figure  7 

Forces  Acting  on  •  Small  Element  of  Cohesionlese  Soil 
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The  arrangement  of  the  stresses  satisfies  the  definition  that  "if  the 
stress  on  a  given  plane  at  a  given  point  is  parallel  to  a  second  plane, 
the  stress  on  the  second  plane  at  the  same  point  must  be  parallel  to 
the  first  plane"  (Ref  50:290).  Rankine  has  shown  that  these  conjugate 
stresses  are  related  to  each  other  by  the  formula 

(7) 

where  F  -  lateral  stress  parallel  to  slope 

P  =  vertical  stress  on  plane  parallel  to  slope 
4  *  angle  of  internal  friction  of  soil 
K  *  conjugate  ratio 

Since  the  vertical  stress  P  is  due  to  tne  weight  of  the  medium  over- 
bur 'Wi,  a  gravity  dependency  for  the  lateral  stress  F  can  be  expressed 
by 

F=pgHK  (8) 

where  p  -  mass  density  of  the  medium 
g  -  acceleration  due  to  gravity 
H  *  depth  of  plane  below  surface 

The  vertical  shearing  resistance  is  a  function  of  the  lateral 
stress  through  the  relation 

V  =  M F  '  <’> 

where  7  >  "ertical  shear  resistance 
fJL  *  coefficient  of  friction 
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Therefore,  ohe  vertical  shear  resistance,  which  must  be  overcome  in 
order  to  accelerate  the  particles,  is  directly  proportional  to  gravity. 
The  constant  force  of  the  expanding  gases  generated  by  the  detonation 
of  identical  explosives  at  the  same  DOB  must  overcome  the  restoring 
force  of  gravity  and  the  vertical  shear  resistance.  For  a  decrease 
in  gravity  there  is  a  decrease  in  the  shear  resistance  on  vertical 
planes  which  allows  the  expanding  gases  to  accelerate  the  particles 
to  higher  velocities  during  the  time  that  the  partic3.es  are  in  contact. 
The  higher  velocities  attained  result  in  higher  initial  velocities  of 
the  particle  when  they  begin  free  flight  as  ejected  particles. 

On  the  basis  of  the  above  explanation  it  is  hypothesized  that 
crater  depths  are  inversely  proportional  to  gravity  when  the  fallback 
due  to  gas  acceleration  is  considered. 

Relative  Importance  of  Processes  as  a  Function  of  DOB 

The  final  size  of  a  terrestrial  crater  is  usually  considered  to 
be  a  function  of  the  type  and  energy  of  the  explosive,  the  medium  in 
which  the  explosion  takes  place,  and  the  depth  of  burst  of  the  explo¬ 
sion.  Since  an  explosive  source,,  constant  in  type  and  energy,  and  a 
homogeneous  cobes ionless  cratering  medium  are  assumed,  only  the  DOB 
is  discussed. 

The  DOB  at  which  the  explosion  takes  place  is  probably  the  most 
important  factor  In  determining  crater  size.  The  reason  for  this  is 
that  the  four  cratering  processes  previously  discussed  are  very  sensi¬ 
tive  to  the  DOB.  This  sensitivity  is  in  turn  reflected  in  any  com¬ 
parison  between  crater  dimensions  and  DOB  such  as  is  illustrated  in 
Fig.  8.  As  is  shown  in  this  figure,  the  depth  and  diameter  of  a 
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peaks  ox  these  representative  curves  are  influenced  by  the  type  of 
medium  in  which  the  crater  is  formed  and  the  explosive  used.  However, 
the  concave  downward  shape  of  these  curves  is  characteristic  of  all 
mediums  and  explosives. 

When  the  charge  is  detonated  on  the  surface  of  the  medium,  com¬ 
paction  and  scouring  are  the  only  processes  available  to  form  the 
crater.  Because  no  portion  of  the  medium  is  located  above  the  charge, 
spalling  and  gas  acceleration  have  no  influence  on  crater  formation. 
For  a  shallow  DOB,  spalling  becomes  the  dominant  process  in  the  forma¬ 
tion  of  a  crater  (Ref  1x1:3395).  Because  of  the  high  velocities  im¬ 
parted  to  the  medium  by  spalling,  the  expanding  gases  escape  into  the 
atmosphere  without  having  to  expend  their  energy  on  the  medium  above 
the  detonated  charge.  Compaction,  and  to  a  lesser  extent,  scouring, 
are  still  involved  in  the  crater  formation.  A r  the  DOB  is  increased, 
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an  optimum  DOB  is  achieved.  At  optimum  DOB,  all  processes  except 
scouring  play  ar.  import  ant,  role-  Gas  acceleration  tends  to  be  the 
dominant  factor,  followed  by  spalling  and  then  compaction.  As  the  DOB 
is  further  increased,  the  relative  importance  of  gas  acceleration  in¬ 
creases.  Spalling  continues  to  decrease  because  the  shock  wave  in¬ 
tensity  is  reduced  by  the  time  it  reaches  the  surface  from  a  deeply 
buried  charge.  The  role  of  compaction  also  continues  to  decrease 
since  the  shear  strength  of  the  medium  increases  with  the  amount  of 
medium  overburden.  The  overall  effect  is  that  the  crater  dimensions 
decrease.  When  the  charge  is  buried  deeply  enough,  the  cratering 
processes  are  so  ineffective  that  all  of  the  disturbed  medium  subsides 
back  into  the  depression,  leaving  no  crater  at  all.  Eventually,  the 
medium  attenuates  the  entire  explosion  to  such  an  extent  that  no  motion 
is  observed  on  the  surface  of  the  medium.  Figure  9,  modified  from 
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Nordyke,  summarizes  the  relative  importance  of  the  cratering  processes 
as  the  EGB  is  increased  (Ref  k2 r 3US? ) •  Based  upon  the  discussion 
contained  in  this  section,  it  can  be  hypothesized  that  the  ordinates 
of  points  on  the  total  dimension  curve  in  Fig.  ?  vary  inversely  with 
gravity. 

Summary 

In  this  section,  it  was  indicated  that  a  change  in  gravity  af¬ 
fects  each  cratering  process.  The  effectiveness  of  the  processes  due 
to  these  changes  and  the  effect  on  the  crater  depth  and  diameter  are 


summarized  below: 

Compaction 

Scouring 

Spalling 

Gas  Acceleration 

Diameter 

Varies 

inversely 

Varies 

inversely 

No 

change 

No 

change 

Depth 

Varies 

inversely 

Varies 

inversely 

Varies 

inversely 

Varies 

inversely 

The  processes  of  compaction  and  scouring  are  most  significant  at 
shallow  depths  of  burst.  At  the  optimum  DOB,  all  cratering  processes 
contribute  to  the  crater  formation.  For  deep  depths  of  burst,  gas 
acceleration  is  the  most  significant  process.  It  is  hypothesized  that 
the  combined  effect  of  the  processes  results  in  crater  dimensions 
varying  inversely  with  gravity. 


III.  Experimental  Procedure 


The  purpose  of  this  experiment  was  to  determine  the  effect  of 
gravity  on  the  formation  of  explosion  craters.  The  experiment  involved 
the  detonation  of  a  small  explosive  in  a  cohesionless  medium  under 
various  gravity  conditions.  This  section  provides  a  brief  explanation 
of  the  manner  in  which  the  experiment  was  conducted.  The  main  steps 
in  this  experiment  were  (l)  securing  the  materials  and  equipment,  (2) 
calibrating  the  equipment,  (3)  conducting  ground  tests,  and  (U)  con¬ 
ducting  flight  tests.  For  a  more  detailed  explanation  of  the  equip¬ 
ment,  calibrations,  sources  of  error,  and  false  leads,  refer  to  Appen¬ 
dix  A. 

Materials  and  Equipment 

The  aircraft  for  varying  gravity  was  available  in  the  Aeronautical 
Systems  Division  at  Wright-Patterson  Air  Force  Base,  Ohio.  The  crater¬ 
ing  medium  selected  was  Flint  Shot  Ottawa  Sand  because  of  its  char¬ 
acteristics  of  uniformity  and  purity  and  its  cohesionless  nature.  Ap¬ 
pendix  B  discusses  in  detail  the  properties  of  the  cratering  medium. 

The  detonations  were  set  off  in  sand  on  the  floor  of  a  plywood  con¬ 
tainer  which  was  U  feet  by  h  feet  on  a  side  and  5  feet  high.  The 
small  explosive  used  was  a  6.0-grain  squib  which  disintegrates  by  a 
deflagration  process.  Photography  was  used  to  record  the  crater 
depths  and  diameters  during  the  experiment. 


Design  and  Calibration 

The  beginning  phases  of  the  experiment  were  concerned  with  the 
calibration  of  the  equipment.  The  purpose  of  the  calibration  was 
(l)  to  insure  that  the  planned  method  of  obtaining  data  was  yielding 
valid  information,  (2)  to  determine  the  degree  of  accuracy  to  which 
measurements  could  be  made,  and  (3)  to  improve  the  efficiency  of  the 
data  gathering  techniques  without  sacrificing  accuracy.  The  factors 
considered  in  calibration  were  (l)  the  area  and  depth  of  sand  neces¬ 
sary  to  minimize  boundary  effects,  (2)  the  best  method  of  placement  of 
the  squib  in  the  sand,  and  (3)  the  best  method  of  extracting  crater 
parameters  from  the  films. 

Ground  Tests 

Ground  tests  were  performed  to  determine  crater  dimensions  under 
1.0  g  conditions.  These  tests  consisted  of  measuring  craters  formed 
under  varying  depths  of  burst.  The  depths  of  burst  were  varied  in 
one-half  inch  intervals  from  surface  bursts  to  a  two  inch  DOB  and  in 
one  inch  intervals  from  a  two  inch  DOB  to  depths  where  no  crater  was 
formed.  Twenty  bursts  were  recorded  at  each  DOB  to  provide  sufficient 
data  for  a  statistical  analysis.  The  steps  in  conducting  a  typical 
crater  forming  cycle  were  (1)  preparing  the  medium,  (2)  setting  the 
charge,  and  (3)  detonating  the  charge  and  recording  the  crater  dimen¬ 
sions. 

Flight  Tests 

The  flight  tests  were  conducted  in  the  same  manner  as  the  ground 
tests  with  the  exception  of  the  changes  in  the  gravity  field.  The 
values  of  gravity  chosen  for  the  experiment  were  0.17,  0.38,  and  2.5  g. 
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These  gravity  levels  could  be  maintained  for  a  period  of  15  seconds. 
The  depths  if  burst  were  varied  in  the  same  manner  as  they  were  during 
the  1.0  g  tests.  Ten  shots  were  taken  for  each  condition  of  grnyi+y 
and  DOB.  The  value  of  gravity  for  each  maneuver  was  recorded  on  an 
oscillograph  for  reference.  In  addition,  five  craters  were  formed  at 
one  and  four  inch  depths  of  burst  in  flight  under  1.0  g  conditions. 
This  was  done  to  determine  any  possible  effect  on  the  crater  dimen¬ 
sions  due  to  aircraft  vibrations. 
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IV.  Discussion  of  Observed  Gravity  Effects 

The  data  from  the  cratering  experiment  is  evaluated  in  this 
section  to  determine  the  effect  of  gravity  on  crater  formation.  The 
experimental  results  are  compared  with  the  previously  suggested  ef¬ 
fects  of  gravity  on  cratering.  Conclusions  are  drawn  concerning  the 
influence  of  gravity  on  crater  formation  for  small  deflagrating  pro¬ 
cesses  in  a  dry  cohesionless  silica  sand. 

Diameters  As  a  Function  of  Depth  of  Burst 

Crater  diameters  can  be  displayed  graphically  as  a  function  of 
the  DOB  of  the  charge.  The  data  on  which  the  graphs  are  based  is 
tabulated  in  Appendix  C.  In  addition,  profiles  of  typical  craters 
formed  at  each  DOB  and  value  of  gravity  are  drawn  in  Appendix  D.  Each 
curve  of  diameter  versus  DOB  is  valid  for  a  given  value  of  gravity. 
Figures  10,  11,  12,  and  13  are  curves  obtained  from  experimental  data 
for  diameters  for  0.17,  O.38,  1.0,  and  2.5  gj  respectively.  Figure 
111  is  a  consolidation  of  these  curves.  Main  points  for  discussion 
concerning  gravity  effects  are  (l)  inherent  scatter  and  possible 
errors,  (C)  significant  differences  in  the  diameters  for  a  given  DOB, 
(3)  changes  in  optimum  DOB,  and  (U)  maximum  DOB  for  which  a  crater  is 
formed. 

In  Section  II,  it  was  hypothesized  that  crater  diameters  vary 
inversely  with  gravity.  Figure  lli  indicates  that  this  hypothesis  may 
be  valid  but  before  the  conclusion  can  be  accepted  the  inherent 
scatter  and  errors  in  the  experimental  data  must  be  checked.  An 
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inspection  of  the  variances  and  ranges  of  the  samples  listed  in 
Table  VI  of  Appendix  C  shows  the  samples  for  O.C,  0.5,  and  5-0  inch 
DOB  at  0.1?  g  have  the  largest  variances.  The  scatter  can  be  at¬ 
tributed  to  aircraft  vibration  and  difficulty  in  measuring  crater 
depths  and  diameters  for  some  craters  of  irregular  shapes. 

Although  vibration  should  affect  shallow  depths  of  burst  at  all 
values  of  gravity,  the  effects  are  more  pronounced  at  the  low  values 
of  gravity.  At  low  values  of  gravity,  the  lateial  pressure  of  the 
sand  and  the  weight  of  the  squib  are  both  decreased.  This  allows 
vibrations  to  move  the  squib  irom  its  desired  position.  During  the 
experiment,  it  was  n^ted  that  the  squio,  when  set,  at  shallow  depths 
of  burst,  had  a  tendency  to  move  upward.  This  was  due  to  tension  in 
the  lead-in  wires  and  firing  wires.  An  effort  was  made  to  eliminate 
this  tension  by  adjustment  of  the  wires,  but  this  was  not  always  ef¬ 
fective.  Squib  movement  from  desired  settings  has  some  influence  on 
crater  diameters  as  indicated  in  Table  VII  of  Appendix  C,  where  di¬ 
ameters  are  compared  when  the  DOB  is  slightly  varied  from  a  given 
position.  Diameters  of  craters  that  are  formed  when  squibs  move  up¬ 
ward  are  smaller  than  diameters  of  crater.,  formed  when  squib  movement 
does  not  occur. 

The  scatter  in  the  data  for  5-0  inch  depths  of  burst  at  0.17  g 
is  caused  by  measurement  difficulties  at  this  condition.  The  final 
crater  shape  was  poorly  defined  and  irregular  because  of  the  deep  DOB 
Craters  were  barely  visible  at  5.0  inch  DOB  at  O.38  g  and  no  craters 
were  formed  at  all  for  this  DOB  at  1.0  and  2.5  g. 

In  summary,  the  charge  aid  not  remain  at  the  desired  DOB  at  shal 
low  depths  of  burst  and  lower  values  of  gravity.  If  the  charge  had 


uhe  ordinate  values  of  the  curves  would 


remained  a:  t:.c  d-s:'red  1.5, 
have  been  son.. ...•*„  higher  than  those  shown  in  Figs.  10  through  lL . 
''rater  shapes  lor  deep  depths  of  burst  were  difficult  to  measure 
which  increased  the  scatter  in  the  data. 

An  understanding  of  the  nature  of  the  scatter  and  errors  in  the 
experiment  permits  a  more  reliable  estimate  to  be  made  concerning 
significant. gravity  effects  on  crater  diameters.  A  statistical  anal¬ 
ysis  of  the  experimental  data  was  made  which  compared  diameters  at 
given  depths  of  burst,  but  which  were  formed  at  two  different  gravity 
values.  Ti.  test  used  is  called  the  two-sample  t-test  and  the  hy¬ 
pothesis  tested  was  that  the  mean  diameters  of  the  compared  samples 
were  equal.  The  alternative  condition  was  that  the  sample  mean  for 
diameters  formed  under  a  given  gravity  field  was  greater  than  the 
sample  mean  for  diameters  formed  under  a  higher  gravity  field.  If 
the  hypothesis  that  the  means  were  equal  could  be  rejected,  gravity 
had  an  effect  crater  formation.  Appendix  E  gives  a  detailed  explan¬ 
ation  of  the  t-test  and  contains  the  tabulated  results  of  this  test. 
This  tabulated  data  reveals  that  the  changes  in  diameters  are  sta¬ 
tistically  significant  for  all  comparisons  except  two. 

One  of  the  comparisons  where  the  hypothesis  could  not  be  rejected 
was  at  0.0  inch  DOB  for  0.17  and  O.38  g.  At  this  DOB,  the  squib  was 
completely  unconfined  by  the  cratering  medium.  The  lack  of  confining 
pressure  on  the  squib  can  result  in  a  large  portion  of  the  charge 
being  scattered  without  burning.  This  property  of  the  squib  is  ex¬ 
plained  in  the  performance  characteristics  section  of  the  specifica¬ 
tions  shown  in  Fig.  4I,  Appendix  A.  The  high  3peed  films  taken  of  the 
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cratering  process  verified  that  portions  of  the  explosive  mixture  were 
scattered  from  the  point  of  detonation  for  surface  bursts.  Such 
scattering  resulted  in  a  reduction  in  the  effective  yield  of  the  ex¬ 
plosive.  The  effective  yield  of  an  explosive  is  defined  here  to  be 
that  portion  of  the  total  energy  yield  which  contributes  to  crater 
formation.  Becauserthe  reduction  in  the  effective  yield  is  an  uncon¬ 
trolled  variable  at  this  DOB,  the  variance  of  the  data  is  increased. 
This  large  variance  tends  to  mask  the  gravity  effects  on  the  cratering 
process  for  0.1?  and  C.38  g.  When  larger  increments  of  gravity  are 
taken  in  the  t-test,  the  gravity  effects  or  crater  dimensions  become 
apparent.  If  an  instantaneous  detonation  could  have  been  achieved  in 
the  experiment,  the  scatter  of  the  explosive  mixture  would  not.  have 
occurred.  One  can  reasonably  speculate  that  the  effective  yield  would 
not  have  been  reduced  and  the  gravity  effects  might  hive  been  sta¬ 
tistically  significant  between  0.1?  and  O.36  g  for  surface  bursts. 

The  other  comparison  for  which  the  hypothesis  cannot  be  rejected 
is  between  0.17  and  O.38  g  at  0.5  inch  DOB.  For  this  DOB,  lack  of 
adequate  confinement  on  the  squib  is  again  an  important  factor.  In 
this  instance,  the  confinement  on  the  squib  is  due  to  the  lateral 
pressure  exerted  on  the  squib  by  the  medium.  According  to  a  duPont 
technical  representative,  the  minimum  condition  for  consistency  of 
the  detonation  process  is  to  have  the  squib  completely  covered  with 
sand  in  a  1.0  g  field  (Ref  lii; .  This  condition  is  achieved  at  0.5 
inch  DOB  for  1.0  or  2.5  g.  However,  for  reduced  gravity,  the  lateral 
pressure  at  a  given  depth  is  reduced.  This  reduction  would  be  espe¬ 
cially  critical  for  0.17  g  and  shallow  depths  of  burst.  The  gravity 
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dependent  confining  pressure  affects  the  explosive  effective  yield  .and 
may  be  responsible  for  the  crossing  in  the  0.17  and  Q.38  g  curves  as 
shown  in  Fig.  lU .  If  an  instantaneous  detonation  could  have  been 
achieved,  the  crossover  in  the  curves  might  not  have  been  present. 

For  all  other  comparisons  in  the  t-test,  statistically  signifi¬ 
cant  changes  were  noted  for  crater  diameters  as  a  function  of  gravity. 
The  experimentally  observed  changes  in  the  diameters  vary  inversely 
with  changes  in  gravity. 

Optimum  and  Maximum  Depth  of  Burst  for  Diameters 

A  point  of  interest  on  Fig.  Hi  is  the  value  for  the  optimum  depth 
of  burst.  The  optimum  depth  of  burst  is  the  DOB  where  the  maximum 
crater  diameter  is  produced  from  the  detonation.  A  change  in  gravity, 
while  changing  crater  diameters,  does  not  appear  to  influence  the  DOB 
for  optimum  diameters.  It  was  hypothesized  in  the  discussion  of  grav¬ 
ity  effects  that  for  constant  depths  of  burst  crater  diameters  are 
influenced  only  by  changes  in  the  compaction  and  scouring  processes. 

If  this  hypothesis  is  true,  the  change  is  shown  in  Fig.  15,  a  modifi¬ 
cation  of  Fig.  9,  A  change  in  the  ordinate  values  cf  the  compaction 
and  scouring  curve  is  reflected  in  a  similar  change  in  the  ordinate 
values  of  the  total  curve.  A  change  of  this  nature  does  not  appear 
to  shift  the  location  of  the  optimum  DOB. 

It  can  be  observed  from  Fig.  II4.  that  the  maximum  DOB  at  which  a 
crater  with  a  measurable  diameter  is  formed  varies  inversely  with 
gravity.  This  observation  may  be  explained  by  the  fact  that  the 
maximum  depth  at  which  shear  failure  can  develop  along  a  plane  varies 
inversely  with  gravity.  This  influences  the  maximum  DOB  at  which  a 


crater  diameter  is  formed. 


Depths  as  a  Function  of  Depth  of  Burst 

The  depth  of  a  crater  can  also  be  analysed  as  a  function  of  the 
DOB  when  gravity  is  varied.  Figures  16>  17,  18,  and  19  are  curves  of 
depth  versus  DOB  for  0.17,  0.38,  1.0,  and  2.£  g,  respectively.  Figure 
20  is  a  consolidation  of  these  curves. 

Some  of  ohe  scatter  in  the  data  for  crater  depths  is  due  to  air¬ 
craft  vibration.  Additionally,  it  is  noted  that  the  curve  in  Fig.  17 
for  O.38  g  has  two  maximum  points.  Of  the  four  curves,  this  is  the 
only  one  which  has  this  characteristic.  A  detailed  investigation  was 
made  to  determine  the  reason  for  this  anomaly.  A  check  of  the  gravity 
values  attained  during  the  flight  and  a  remeasurement  of  the  craters 
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inversely  with  gravity,  the  maximum  DOB  also  varies  inversely  with 

gravity. 

Relationship  Between  Dimensions  and  Gravity 

At  this  point,  it  has  been  established  that  crater  dimensions 
vaiy  Inversely  with  gravity.  The  next  step  is  to  try  to  determine  a 
mathematical  relationship  for  the  dimensions  as  a  function  of  gravity. 

The  dimensions  should  have  some  finite  value  even  as  the  gra  ufcy 
field  approaches  zero.  The  dimensions  of  craters  formed  at  very  low 
values  of  g  are  determined  by  the  manner  in  which  the  sand  dissipates 
the  energy  of  the  detonation  through  inertial  resistance. 

Examining  the  limiting  condition  when  the  gravity  field  becomes 
very  large,  one  expects  the  crater  dimensions  to  approach  zero  a*'  a 


1*6 


limit.  This  is  cased  upc:.  the  inability  of  the  energy  '  f  the  explo¬ 
sion  to  overcome  the  large  shear  resistance  which  i3  present.  In 
addition.,  the  range  of  any  particle  which  is  given  an  initial  velocity 
becomes  negligible  as  is  illustrated  by  Eq.  (6).  Considering  the  above 
limiting  cases  whiol.  must  be  satisfied.,  an  exponential  function  for 
the  dimension  in  terms  of  gravity  appears  to  be  a  suitable  relationship 
to  assume. 

Plots  of  diameters  and  depths  as  a  function  of  gravity  for  a 
given  DOB  are  shown  in  Figs.  22  through  35.  It  can  be  seen  that  the 
sample  points  do  indeed  exhibit  a  trend  toward  a  a'icreasing  exponen¬ 
tial  function  as  the  gravity  field  is  increased.  In  order  to  obtain 
parameters  so  that  gravity  effects  could  be  analyzed,  an  exponential 
function  was  assumed  of  the  form 

D  -  Ae~gB  (10) 

where  D  -  dimension  of  the  crater 
A  =  arbitrary  constant 

B  *  measure  of  sensitivity  of  dimension  to  gravity  changes 
g  *  acceleration  due  to  gravity 

Thii  parameters  A  and  Z  .  aJ;u3t»a  • .  achieve  the  best  to  the 
experimental  data  in  the  least-squares  sense.  The  effectiveness  of 
this  curve  fitting  technique  was  limited  by  the  availability  of  only 
four  data  points  for  each  curve.  Values  of  A  and  B  obtained  in  the 
curve  fitting  calculations  and  the  least-squares  deviations  are  listed 
in  Table  I.  In  addition,  Figs.  22  through  35  show  plots  of  Eq.  (10) 
for  the  different  depths  of  burst  using  values  of  A  nrd  B  as  calculated. 
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TABLE  I 

VALUES  OF  A  AND  B  AS  A  FUNCTION  OF 
DEPTH  OF  BURST  AND  CRATER  DIMENSIONS 


DOB 

in 

inches 


0.0 

0.5 

1.0 

1.5 

2.0 

3.0 

4.0 


B 


A 


Diameter 


Least-Squares 
Deviation  * 


0.116 

24 

0.121 

32 

0.144 

36 

0.150 

38 

0.151 

37 

0.160 

33 

0.168 

28 

•4529  0.1581 
•8782  2.0721 
•9874  0.0844 
•3664  0.4196 
•6850  0.2413 
•3658  0.0924 
.1812  0.0582 


Depth 


0.0 

0.075 

0.5 

0.042 

1.0 

0.085 

1.5 

0.156 

2.0 

0.209 

3.0 

0.247 

4.0 

0.225 

4.4240  0.0339 
5*8720  0.2882 
6.5554  0.0087 
6.3021  0.0309 
5.2473  0.0074 
3.5691  0.0014 
2.1196  0.0075 


*  The  deviation  is  the  value  of  the  sua  of  the  squares  of  the 
differences  between  the  set  of  data  points  for  the  dimensions 
at  a  given  DOB  and  the  assumed  exponential  curve,  D  -  Ae**P. 
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The  magnitude  of  the  deviation  from  a  perfect  fit  is  the  determining 
factor  in  the  least-squares  technique.  Table  I  and  the  graphs  show 
.hat  all  the  deviations  are  less  than  1.0  except  at  0.5  inches  DOB  for 
crater  diameters.  The  dependence  of  effective  yields  on  gravity 
causes  this  large  deviation.  In  addition  values  of  B  change  with  depth 
of  burst.  These  trends  are  illustrated  in  Figs.  36  and  37,  which  are 
consolidations  of  the  plots  of  Eq.  (10).  In  these  figures,  the  value 
of  B  is  reflected  in  the  slopes  of  the  curves. 

The  range  of  B  is  small  for  diameters  when  compared  to  the  range 
of  B  for  crater  depths.  This  is  expected  because  only  compaction  and 
scouring  cause  incremental  changes  in  the  diameter  when  gravity  is 
varied.  The  B  values  slowly  increase  as  the  DOB  increases  from  0.C  to 
U.O  inches.  The  increase  means  that  the  relative  importance  of  com¬ 
paction  and  scouring  due  to  changes  in  gravity  increases  with  in¬ 
creasing  DOB. 

The  larger  range  of  B  for  depths  than  for  diameters  may  be  ex¬ 
plained  by  the  hypothesis  that  all  cratering  processes  cause  changes 
in  the  depth  of  a  drater  when  gravity  is  varied.  The  relatively  low 
values  of  B  at  0.0,  0.5,  «nd  1*0  inch  DOB  are  caused  rf  the  low  angle 
of  ejection  and  high  initial  velocities  of  the  ejected  particles  which 
throw  them  clear  of  the  crater  regardless  of  the  value  of  gravity.  At 
depths  of  burst  less  then  1.0  inches,  the  change  in  depth  is  determined 
mostly  by  compaction.  The  particularly  low  value  of  B  at  0.5  inch 
DOB  reflects  the  influence  of  gravity  on  the  effective  yield  of  the 
explosive.  The  increase  in  the  fallback  with  increasing  DOB  is  evi¬ 
denced  by  larger  values  of  B  at  1.5,  2.0,  and  3*0  inch  DOB.  The 
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large  increase  in  B  with  greater  depths  of  burst  tends  to  support  the 
hypothesis  that  the  g as  acceleration  process  is  more  influenced  by 
gravity  changes  than  the  other  processes  when  depths  are  considered. 
This  hypothesis  was  put  forth  as  an  explanation  for  the  shift  in  the 
optimum  DOB  for  depths.  The  B  values  reach  a  maximum  and  decrease  in 
the  vicinity  of  3»0  inch  DOB.  This  decrease  occurs  because  the  angle 
of  ejection  of  particles  at  depths  of  burst  greater  than  3*0  inches 
become  so  nearly  vertical  that  the  effect  of  gravity  on  particle  ranges 
is  reduced. 

It  should  be  emphasized  here  that  the  exponential  expression 
which  was  chosen  to  represent  crater  dimensions  as  a  function  of  grav¬ 
ity  may  not  be  the  correct  expression  since  a  complete  theoretical 
analysis  was  not  developed  which  proves  the  validity  of  this  relation¬ 
ship.  However }  this  exponential  expression  approximates  the  data  well 
enough  to  determine  that  the  sensitivity  of  crater  dimensions  to 
changes  in  gravity  is  a  function  of  the  DOB. 
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-  V.  Conclusions  and  Recommendations 

Conclusions 

The  main  points  to  be  made  from  this  experiment  are: 

1.  Crater  parameters  which  vary  inversely  with  gravity  are: 

a.  crater  depths  and.  diameters 

b.  optimum  depth  of  burs*-' for  crater  depths 

c.  maximum  depth  of  burst  for  which  a  visible  crater  is 

formed. 

2.  The  optimum  depth  of  burst  for  crater  diameters  is  not  a 
function  of  gravity. 

3.  The  sensitivity  of  changes  in  crater  dimensions  to  changes 
in  gravity  is  a  function  of  the  depth  of  burst. 

Recommendations 

The  following  topics  in  the  field  of  cratering  in  varied  gravity 
are  recommended  for  further  study: 

1.  Use  of  an  explosive  with  a  higher  detonating  velocity  to 
compare  results  from  a  deflagrating  explosive. 

2.  High  speed  analysis  of  the  cratering  process  to  include 
calculations  of  throvout  patterns. 

3.  Cratering  in  different  mediums. 

U.  Measurement  of  crater  volumes. 

5.  Quantitative  theoretical  study  of  cratering  in  varied  gravity 
fields  based  on  the  discussion  developed  in  this  paper. 
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Appendix  A 

Details  of  Experimental  Procedure 
Equipment  and  Materials 

The  equipment  and  materials  required  for  the  experiment  were  a 
special  aircraft  for  varying  gravity  conditions,  a  cohesionless  sand, 
a  container  in  which  the  cratering  process  could  take  place,  a  small 
explosive  source,  and  devices  to  measure  the  crater  parameters. 

The  Zero-Gravity  Aircraft.  The  special  aircraft,  a  twin-engined, 
propeller-driven.  Air  Force  C-131B,  was  available  in  the  Aeronautical 
Systems  Division  at  Wright-Patterson  Air  Force  Base,  Ohio.  Gravity 
was  varied  under  controlled  conditions  during  the  flights.  Figure  38 
shows  the  nose  of  the  aircraft  and  a  flight  profile  which  indicates 
the  details  of  the  maneuver  performed  to  vary  gravity  inside  the  air¬ 
craft.  By  flying  an  electronically  programmed  maneuver  similar  to  that 
for  the  aero  gravity  condition  shown  in  Fig.  38,  any  value  of  g  up  to 
1.0  g  could  be  attained  for  periods  up  to  15  seconds.  Values  greater 
than  1.0  g  and  up  to  2.5  g  were  achieved  by  performing  tight  tum;i 
instead  of  the  parabolic  maneuver.  Details  about  the  aircraft  and  its 
flight  plan  can  be  found  in  Reference  3* 

The  Cratering  Medium.  Because  of  its  homogeneous  qualities  com¬ 
pared  to  other  soils,  e  dry,  cohesionless,  silica  sand  was  selected  as 
the  medium  in  which  to  fora  the  craters.  Flint  Shot  Ottawa  Sand  was 
the  brand  picked  because  of  its  availability  and  outstanding  character¬ 
istics  of  uniformity  and  purity.  Figure  39  illustrates  these 
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characteristics  by  means  of  a  close-up  photograph  of  the  sand  particles. 
Based  upon  a  number  of  soil  mechanics  tests,  it  was  determined  that  the 
sand  maintained  consistent  values  of  density  and  moisture  content 
throughout  the  experiment.  A  discussion  of  these  tests  and  their  re¬ 
sults  are  described  in  Appendix  B. 

The  Sand  Container.  A  container  was  necessary  in  order  to  confine 
the  cratering  process  during  the  experiment.  This  was  especially  im¬ 
portant  on  the  aircraft  as  sand  could  not  be  permitted  to  come  in  con¬ 
tact  with  the  mechanical  and  electrical  equipment  in  the  aircraft. 

The  box  was  designed  to  be  large  enough  so  that  it  would  not  influence 
the  cratering  process  and  small  enough  to  move  it  on  and  off  the  plane. 
Considering  these  criteria,  the  box  was  constructed  from  3/U-inc.h  ply¬ 
wood  with  a  h  feet  by  U  feet  base.  The  height  of  the  box  was  5  feet. 

One  wall  of  the  box  was  constructed  as  a  door  to  allow  access  inside 
the  box.  In  addition,  a  plexiglass  section  was  mounted  on  tcp  to  al¬ 
low  observation  of  the  cratering  process.  Figure  UO  is  a  photograph 
of  the  container. 

The  Small  Explosive  Source.  The  dimensions  of  the  sand  container 
were  primarily  determined  by  the  interior  dimensions  of  the  aircraft, 
this  placed  an  upper  limit  on  the  magnitude  of  the  explosive  that, 
could  be  used  for  the  cratering  process .  An  explosion  of  too  great  a 
magnitude  would  cause  particles  to  ro bound  from  the  walls  of  the  con¬ 
tainer  end  influence  the  final'  crater  form  by  falling  back  into  the 
crater  when  they  would  normally  be  throw  clear.  The  possibility  of  a 
reflected  shock  wave  from  the  walls  of  the  container  Interacting  with 
the  normal  cratering  process  also  had  to  be  avoided.  In  addition. 


Figure  UC 
Sand  Container 


flight  safety  demanded  a  detonation  of  small  magnitude.  A  6.0  grain 
explosive  squib  was  found  to  be  suitable  as  the  explosive  device  under 
the  described  restrictions.  Use  squibs  that  were  used  are  electrically 
detonated  explosives  normally  used  to  actuate  release  mechanisms  such 
as  wing  tip  tank  releases  on  Air  Force  aircraft.  Figures  1*1  and  li2 
show  the  details  of  the  squib  used  in  the  experiment.  The  load  toler¬ 
ance  of  the  squib  is  6.0  -  0.U  grains  of  explosive.  The  duPont 
50/2^/25  mixture  in  the  3quibs  does  not  detonate  instantaneously  as  is 
the  case  with  an  explosive  like  TNT.  Instead,  the  mixture  burns  rapid¬ 
ly  to  produce  high  pressure  levels.  This  type  of  mixture  is  known  as 
a  deflagrating  compound  and  is  in  a  different  class  from  TNT.  TNT 
expends  its  energy  with  a  strong  shock  wave  whereas  the  duPont  50/25/25 
mixture  expends  its  energy  through  a  relatively  slow  build-up  of 
pressure  over  a  longer  period  of  time  than  TNT.  This  results  in  a 
weaker  shock  wave. 

Photographic  Measurements.  Photographic  means  were  used  to 
measure  the  era  ter  parameters  during  the  experiment.  The  crater  par¬ 
ameters  chosen  for  observation  were  crater  depth  and  diameter.  These 
dimensions  are  indicated  on  Fig.  1.  Photography  was  chosen  because 
aircraft  vibrations  disturbed  craters  before  manual  measurements  could 
be  made.  This  also  allowed  the  maximum  number  of  detonations  to  be 
performed  dnidng  the  limited  flight  time  available  for  the  experimental 
phase.  A  16-ram  motion  picture  camera  of  the  B1A  type  was  used  with  a 
13-mm  lens  and  a  frame  rate  of  6k  frames  per  second.  The  camera  was 
electrically  powered  with  a  28- volt  direct  current  motor.  Figure  U3 
shows  the  camera  as  it  was  mounted  on  the  sand  container  in  a  position 
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to  record  the  cratering  process.  The  electronic  equipment  that  con¬ 
trolled  the  camera  and  cratering  process  is  also  shown  in  this  figure. 
This  circuitry  enabled  the  experiment  to  be  controlled  by  the  activa¬ 
tion  of  a  single  push-button  switch.  Pushing  this  button  started  the 
camera  and  flood  light.  After  a  2  second  delay,  28- volts  were  auto¬ 
matically  sent  to  the  prepared  squib.  After  initiation  of  the  squib, 
the  camera  and  light  continued  to  function  for  3  seconds.  The  filming 
and  lighting  sequence  then  automatically  terminated  and  a  new  crater¬ 
ing  sequence  could  be  started.  A  wiring  diagram  for  the  circuit  is 
3hown  in  Fig.  hh-  The  lighting  apparatus  is  shown  in  Fig.  U5.  These 
lights  were  energized  by  110- volt  alternating  current. 

In  addition,  high-speed  photography  was  used  to  record  the  motion 
of  the  sand  during  the  cratering  process.  A  high  speed  sequence  was 
taken  of  one  crater  for  each  condition  of  gravity  and  DOB.  This  was 
done  to  provide  an  additional  means  of  analyzing  the  processes  occur¬ 
ring  during  crater  formations.  The  camera  used  was  a  16-m  Fairchild 
Motion  Analysis  Camera,  Model  Number  MS  101,  powered  by  a  28- volt, 
direct  current  motor.  A  13-nm  lens  was  provided  on  the  camera.  The 
camera,  which  operated  at  a  maximum  rate  of  2000  frames  per  second, 
is  shown  in  Fig.  1*6. 

Calibration  Teats 

The  calibration  phase  of  the  experiment  in  this  thesis  included 
determinations  of  (1)  the  boundary  effects  due  to  the  container  and 
depth  wire,  (2)  a  satisfactory  method  of  placement  of  the  squib  in  the 
•and,  and  (3)  satisfactory  methods  of  extracting  crater  parameters 
from  the  films. 
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Boundary  Effects.  In  order  to  insure  that  the  cratering  process 
was  not  influenced  by  any  reflections  of  the  shock  wave  from  the 
..ottom  of  the  container,  the  depth  of  sand  used  as  the  medium  was 
carefully  investigated.  A  comparison  of  crater  parameters  in  dif¬ 
ferent  depths  of  sand  indicated  whether  or  not  boundary  effects  were 
influencing  crater  dimensions.  It  was  felt  that  enough  sand  between 
the  explosion  and  the  container  would  damp  out  any  reflected  shock 
waves  that  would  adversely  affect  the  crater.  The  different  depths 
investigated  were  5,  10,  12,  13,  and  1U  inches  of  sand.  Different 
size  craters  were  formed  between  5,  10,  12,  and  13  inches  of  sand  for 
a  constant  source  of  energy  and  DOB.  The  differences  were  felt  to  be 
due  to  boundary  effects  from  the  container.  A  statistical  analysis 
using  the  two-sample  t-test  indicated  that  there  was  not  a  significant 
difference  in  compared  dimensions  of  craters  formed  in  13  end  1U  inches 
of  sand  (see  Appendix  E) .  Therefore  negligible  boundary  effects  were 
present  with  lU  inches  of  sand  in  t!  s  container.  Since  there  was 
more  than  lh  inches  of  sand  between  the  point  of  detonation  and  the 
sides  of  the  container,  the  assumption  was  made  that  there  were  no 
boundary  effects  from  the  side  of  the  container. 

The  thin  wire  shown  in  Pig.  U5  was  used  as  an  aid  in  extracting 
the  crater  dimensions  from  the  motion  picture  film.  This  technique 
will  be  explained  in  the  discussion  of  crater  measurements  contained 
in  this  section.  The  wire  did  not  influence  the  crater  dimensions  as 
comparisons  of  depths  and  diameters  of  craters  formed  under  the  same 
conditions  except  with  or  without  the  wire  showed  no  statistically 
significant  difference. 
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Squib  Placement.  Several  methods  of  placement  of  the  squib  to 
obtain  the  desired  DOB  were  tested.  The  squib  had  to  be  (1)  emplaced 
quickly  to  conserve  time  on  the  aircraft,  (2)  attached  securely  to 
withstand  aircraft  vibrations,  and  (3)  positioned  accurately  at  a  given 
depth.  It  was  determined  that  a  plarvirr,  >f  the  squib  with  the  wires 
pointing  upward  worked  most  effectively.  The  center  of  detonation  was 
estimated  from  Fig.  Ul.  Based  upon  this  figure,  the  midpoint  distance 
between  the  end  of  the  squib  and  the  first  crimp  was  chosen  as  the 
center  of  the  detonation.  Distances  were  measured  upward  from  this 
point  and  the  desired  depth  settings  were  marked  on  the  wires.  After 
the  sand  was  leveled  to  a  depth  of  li*  inches,  the  squib  was  inserted 
into  the  sand  at  a  point  next  to  the  depth  wire  with  the  lead-in  wires 
pointing  upward.  In  order  to  insure  that  the  squib  was  vertical  and 
that  the  lead-in  wires  were  straight,  the  squib  was  inserted  to  a  depth 
below  the  final  desired  depth  and  then  pulled  upward  to  the  final  po¬ 
sition.  The  final  position  was  noted  by  means  of  a  mark  painted  on  the 
lead-in  wires.  The  above  procedure  proved  to  be  the  most  accurate 
method  of  achieving  the  correct  DOB.  After  the  squib  was  attached  tc 
the  firing  wires,  the  container  was  closed  and  the  firing  sequence  was 
initiated. 

For  depths  of  burst  greater  than  two  inches,  it  was  necessary  to 
slightly  modify  the  above  procedure  due  to  the  inability  of  inserting 
the  squib  to  the  desired  depth  without  bending  the  lead-in  vires.  For 
these  depths,  a  depression  was  made  in  the  medium  where  the  squib  was 
to  be  emplaced.  This  technique  is  shown  in  Fig.  U7.  The  depression 
allowed  the  squib  to  be  placed  at  the  desired  depth  without  having  to 
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be  forced  through  more  than  two  inches  of  the  medium.  The  depression 
was  then  filled  in  and  the  sand  was  leveled.  Final  positioning  of  the 
squib  was  accomplished  by  pulling  the  lead-in  vires  upward  to  insure 
the  proper  DOB  and  vertical  attitude.  Experimentation  with  this  method 
indicated  that  the  squib  remained  close  to  vertical  and  that  the  depth 
could  be  consistently  set  to  within  1/6  inch  of  the  desired  depth. 

Some  error  in  the  data  could  have  resulted  from  this  method  of  3quib 
placement ,  but  the  error  was  random  rather  than  systematic. 

Crater  Dimension  Measurements .  The  dimensions  of  the  craters 
formed  during  this  experiment  were  obtained  from  motion  picture  film. 

The  validity  and  accuracy  of  the  measurements  taken  from  the  film  were 
checked  by  comparing  the  film  measurements  with  manual  measurements 
taken  of  the  same  crater. 

The  crater  measurements  were  extracted  from  the  film  with  the  aid 
of  a  thin  wire  and  a  photographed  scale.  The  wire  :7a 3  placed  in  the 
container  as  shown  in  Fig.  U8  and  a  reference  mark  was  painted  on  the 
wire.  By  observing  on  the  film  the  intersection  of  this  wire  with  -the 
sand  before  and  after  the  detonation,  the  depth  of  the  crater  was  de¬ 
termined.  The  diameter  was  directly  observed  on  the  film  as  the  rim- 
to-rira  distance  normal  to  the  line  of  sight  of  the  photograph.  A 
photographed  scale  was  required  to  convert  the  scaled  depth  and  diameter 
as  seen  on  the  film  to  true  depth  and  diameter.  This  was  done  with 
the  scale  that  is  shown  in  Fig.  1*9.  This  scale  was  formed  by  placing 
a  vertical  ruler  along  the  depth  wire  with  the  12-inch  line  adjacent 
to  the  top  of  the  reference  mark  on  the  wire.  The  horizontal  ruler 
was  placed  at  the  oilginal  ground  level.  Several  feet  of  film  were 
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Pre-Detonation  Condition 


Figure  k9 

Scale  Used  for  Crater  Measurements 

taken  of  this  scale  from  the  same  location  and  with  the  same  camera  as 
used  to  photograph  the  craters. 

The  actual  crater  dimensions  were  obtained  from  the  film  by  pro¬ 
jecting  the  cratering  sequence  on  a  poster  board  screen^  The  pre- 
detonation  conditions  were  recorded  by  marking  the  position  of  the 
wire  reference  marie  and  the  original  ground  level  on  the  screen.  The 
Intersection  of  the  wire  and  the  sand  clearly  located  the  original 
ground  level.  The  film  was  advanced  to  the  post- detonation  position 
where  the  final  crater  was  recorded.  After  relocating  and  super¬ 
imposing  the  reference  marks  on  the  screen  and  film,  the  wire-sand 
intersection  could  again  be  marked.  The  wire-sand  intersection  is 
seen  in  Fig.  J> 0 .  The  distance  between  this  mark  and  the  original 
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Figure  50 

Final  Crater  and  Wire-Sand  Intersection 


53 


ground  level  was  the  scaled  crater  depth.  The  normal- to- the-line-of- 
sight  rim  separation  distance  across  the  center  of  the  crater  was 
marked  as  the  scaled  diameter.  Figure  1  shows  these  dimensions.  By 
projecting  the  photographed  scale  on  the  screen  and  aligning  the 
reference  marks  again,  the  actual  depths  and  diameters  could  be  read 
directly.  This  procedure  was  further  expedited  by  tracing  the  pro¬ 
jected  scale  on  acetate.  The  acetate  was  then  placed  over  the  screen 
and  properly  aligned  to  measure  the  actual  crater  dimensions.  This 
system  was  valid  as  long  as  the  screen-projector  distance  remained 
constant  for  the  entire  measuring  procedure.  This  fixed  distance  was 
carefully  maintained  during  the  measurement  of  the  craters. 

The  validity  of  this  method  of  measurement  was  proven  as  the 
manual  and  film  measurements  agreed  on  the  average  to  within  -  1/8 
inch  for  the  depths  and  -  lA  inch  for  the  diameter.  In  addition, 
some  of  the  craters  were  remeasured  from  the  film  at  a  later  date  as 
a  check  of  the  consistency  of  the  measuring  technique.  The  average 
change  in  the  crater  dimensions  based  upon  the  remeasured  dimensions 
was  -  lA  inch  for  the  diameter  and  -  1/8  inch  for  the  depth.  The 
limit  of  measurement  of  the  photographed  scale  could  be  interpolated 
to  the  nearest  1/8  inch  when  projected  on  the  screen. 

Sources  of  Error 

There  were  many  possible  sources  of  error  in  the  experiment.  The 
load  tolerance  of  the  squib  was  6.0  i  0.4  grains,  which  allowed  an 
approximate  1%  deviation  in  the  explosive  yield.  The  squib  was  con¬ 
structed  so  that  the  detonation  was  not  spherically  symmetrical.  This 
placed  additional  stress  on  the  necessity  for  accurate  placement 
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techniques .  The  small  magnitude  of  the  yield  caused  a  small  error  in 
placement  to  be  greatly  magnified.  In  Appendix  C  are  tabulated  the 
dimensions  of  craters  formed  at  a  DOB  of  slightly  less  than  0.5  inches 
and  slightly  greater  than  0.5  inches.  The  tabulated  dimensions  show 
that  even  a  small  error  in  placement  can  affect  the  crater  dimensions. 

An  additional  factor  with  the  squib  was  the  nature  of  its  detona¬ 
tion.  Since  the  squib  oumed  rapidly  instead  of  decomposing  instan¬ 
taneously,  the  lack  of  confinement  of  the  charge  at  surface  bursts  per¬ 
mitted  some  scattering  of  unbumed  or  burning  compound. 

Aircraft  vibrations  had  some  influence  on  the  location  of  the 
squib  between  the  time  that  it  was  placed  and  detonated.  Vibrations 
may  have  also  had  some  effect  on  crater  formation.  A  t-test  comparison 
of  crater  dimensions  for  1.0  g  conditions  on  and  off  the  aircraft  is 
shown  in  the  tabulated  data  in  Appendix  £.  The  results  of  this  test 
indicate  that  the  diameters  are  only  slightly  influenced  by  vibrations 
but  the  depths  show  significant  vibration  changes.  The  error  due  to 
vibrations  increase  with  increasing  DOB. 

The  value  of  gravity  in  the  aircraft  was  subject  to  variations  of 
±  .01  g  for  any  particular  maneuver  due  to  turbulence  and  limitations 
of  pilot  control.  It  was  observed  during  the  flights  that  occasional 
longitudinal  acceleration  of  the  aircraft  was  present  during  some  of 
the  maneuvers.  This  was  due  to  Improper  throttle  manipulations  and, 
when  noticed  by  the  crew  or  technicians,  was  corrected.  The  presence 
of  longitudinal  acceleration  primarily  affected  the  depth  of  the 
crater  rather  than  the  disaster.  At  shallow  depths  of  burst,  more  of 
the  medium  would  slide  back  into  the  center  of  the  crater  after  the 
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detonation  process  was  completed.  For  deep  depths  of  burst,  longi¬ 
tudinal  acceleration  of  the  aircraft  resulted  in  a  shift  of  the  fall¬ 
back  pattern.  This  resulted  in  a  shift  in  the  location  of  the  crater 
center  away  from  the  depth  wire.  It  was  felt  that  longitudinal  accel¬ 
eration  was  the  largest  source  of  error  in  the  flight  tests.  In 
addition,  there  were  random  errors  in  measurements  due  to  the  limits 
of  resolution  of  the  film. 

False  Leads 

False  leads  encountered  during  the  experiment  are  listed  here 
briefly  in  hopes  that  they  might  possibly  help  future  investigators 
save  some  time  and  effort  in  -heir  research. 

Other  explosive  sources  than  squibs  were  investigated.  Exploding 
wires  were  considered  but  not  used  as  the  cost  of  constructing  a  de¬ 
vice  to  explode  wires  was  prohibitive  when  compared  to  squibs.  Ad¬ 
vantages  of  exploding  wires  are  ease  of  control  and  safety.  The  ex¬ 
plosion  of  a  wire  is  not  spherically  symmetrical,  however,  and  the 
energy  partition  is  different  from  a  chemical  explosion. 

Crater  parameters  other  than  depth  and  diameter  were  considered. 
An  attempt  to  ascertain  crater  volumes  proved  fruitless  from  the  films 
as  the  diameter  of  the  crater  was  measured  from  the  top  of  the  crater 
rim  rather  than  along  the  original  ground  level.  The  top  of  the 
crater  rim  was  the  only  observable  position  from  which  to  make  a 
diameter  determination  from  the  film.  The  error  in  the  measured 
volume  when  the  rim- to- rim  diameter  was  used  instead  of  the  proper 
ground  level  diameter  ranged  as  high  as  lOOjf.  Rim  heights  could  not 
be  measured  because  of  the  lack  of  reliable  references  on  the  films 
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and  the  small  rim  height.  Rim  heights  seldom  ranged  over  1/2  inch. 

An  attempt  to  make  depth  determinations  by  painting  the  depth 
wire  in  known  increments  was  made.  By  painting  the  wire  in  known  in¬ 
crements,  the  depth  of  the  crater  could  have  been  determined  by  simply 
counting  the  difference  in  the  number  of  increments  visible  before  and 
after  detonation.  The  paint  could  not  withstand  the  explosive  force 
and  the  increments  were  either  obscured  or  destroyed.  The  photographic 
scale  was  devised  as  an  alternate  means  of  measuring  the  crater  depths. 
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Appendix  B 

Properties  of  the  Cratering  Medium 

Ottawa  sand  is  a  high  grade  white  silica  sand  which  is  mined  ex¬ 
clusively  from  a  silica  sandstone  of  exceptional  purity  located  in 
Ottawa,  Illinois.  Because  of  the  high  degree  of  uniformity  and  clean¬ 
liness,  it  has  found  extensive  use  as  a  standard  for  testing.  In  par¬ 
ticular,  the  United  States  Government  uses  certain  grades  of  Ottawa 
sand  in  calculating  the  energy  content  of  chemical  explosives 
(Ref  U3*59) . 

Because  the  properties  of  the  medium  in  which  the  crater  is 
formed  influence  the  size  of  the  crater,  it  was  imperative  that  these 
properties  be  known  and  that  the  variance  in  certain  critical  char¬ 
acteristics  such  as  density  and  moisture  content  be  limited  throughout 
the  experiment.  In  order  to  obtain  this  information,  certain  soil 
tests  were  performed.  The  results  of  these  tests  are  discussed  in  the 
following  paragraphs. 

Grain  Size  Analysis 

The  process  of  separating  a  soil  into  particle-size  groups  by 
shaking  it  through  a  stack  of  wire  screens  with  openings  of  known 
sizes  is  called  sieve  analysis.  The  definition  of  particle  size  is 
the  size  of  a  square  opening  on  which  the  grain  of  soil  is  retained. 
The  recommended  procedure  which  is  outlined  in  Section  IV  of  Reference 
31  was  followed  for  the  conduct  of  the  test. 


98 


Mech/QSF-65-35 


The  result  of  the  test  is  represented  graphically  in  Fig.  $1  by 
means  of  a  grain  size  distribution  curve.  From  the  figure  it  is  evi¬ 
dent  that  the  particles  of  the  soil  sample  used  as  the  cratering 
medium  are  uniform  in  size.  A  further  indication  of  the  general 
nature  of  the  medium  is  given  by  the  terns  known  as  the  coefficient  of 
gradation  and  the  uniformity  coefficient.  The  coefficient  of  grada¬ 
tion  is  defined  as 


while  the  uniformity  coefficient  is  defined  as 

Cu  ■  d60  (12) 

where  the  Dj_  refers  to  the  grain  size  which  corresponds  to  the  i$ 
finer  by  weight  point  on  a  grain  size  distribution  curve,  in  order  to 
be  classified  as  being  uniformly  graded,  the  value  of  the  uniformity 
coefficient  must  be  less  than  U  while  the  value  of  the  coefficient  of 
gradation  must  be  less  than  1  or  greater  than  3  (Ref  56*281).  For  the 
Ottawa  sand  sample,  the  values  of  the  two  coefficients  are  Cg  *  .92 
and  -  1.1*7.  Baaed  upon  these  coefficients,  the  msdium  1s  uniformly 
graded.  For  purposes  of  the  investigation,  a  uniform  soil  sample  was 
desired  since  this  should  decrease  the  variance  in  the  density  of  the 
soil. 


Density 

Tests  were  conducted  in  order  to  determine  if  the  density  of  the 
sample  would  exhibit  only  a  limited  variance.  The  first  test  consisted 
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of  placing  the  dry  Ottawa  sand  in  a  standard  compaction  mold  normally 
used  for  the  Proctor  compaction  test.  The  sand  was  placed  in  as  loose 
a  manner  as  possible.  Using  '■he  equation 

7d  -  «  <«> 

"ra 

where  W  *  weight  of  sand  in  mold  in  pounds 

Vm  *  volume  of  the  mold  which  was  1/30  of  a  cubic  foot 
^  *  density  of  the  sample  in  pounds  per  cubic  foot  (pcf) 
a  density  of  96.8  pcf  was  attained. 

The  second  tes ;  consisted  of  placing  the  dry  Ottawa  sand  in  the 
mold  in  a  loose  state  and  than  subjecting  the  mold  to  a  small  amount 
of  agitation  and  vibration  such  as  it  would  be  subjected  to  during  the 
ground  and  flight  testa.  For  this  test,  the  density  was  calculated  to 
be  10?. 5  pcf. 

The  final  test  consisted  of  performing  a  modified  Proctor  test  on 
the  dry  sample  in  which  the  sand  was  placed  In  the  mold  in  five  layers 
with  each  layer  being  compacted  by  means  of  2$  blows  with  a  ten  pound 
hammer  falling  through  a  distance  of  18  inches.  A  detailed  explana¬ 
tion  of  such  a  test  can  be  found  in  Chapter  7  of  Reference  31.  For 
this  test  the  density  was  calculated  to  be  108  pcf. 

Baseu  upon  the  above  tests,  it  was  concluded  that  the  variation 
in  density  due  to  the  conditions  of  the  experiment  would  be  negligible. 

Moisture  Content 

Since  the  assumption  of  negligible  variation  in  density  depended 
upon  a  negligible  variation  in  the  moisture  content  of  the  sand,  tests 


ditions , 


In  the  first  test,  the  moisture  content  of  sand  which  had  been 
stored  in  its  original  paper  sack  was  determined  according  to  the  pro¬ 
cedures  outlined  in  Chapter  III  of  Reference  31.  The  results  of  the 
tests  on  five  different  samples  indicated  zero  moisture  content. 

A  second  moisture  content  determination  was  made  on  five  samples 
of  sand  which  were  being  used  for  the  ground  tests .  In  this  case,  the 
medium  was  stored  in  a  heated  building.  Again  zero  moisture  content 
was  recorded. 

A  final  test  was  made  on  five  samples  of  sand  which  were  stored 
in  open  containers  in  an  unheated  garage  during  three  days  of  inter¬ 
mittent  rain.  During  this  time,  the  door  of  the  garage  remained  open 
2h  hours  a  day  in  order  to  make  the  air  in  the  garage  as  humid  as 
possible.  The  moisture  content  for  three  of  the  samples  was  rerc. 

For  the  other  two  samples,  the  moisture  content  was  0.l£.  No  water 
was  ever  added  to  the  sand  during  the  moisture  content  tests  since  the 
sand  and  the  box  in  which  it  was  contained  were  never  exposed  to  any 
moisture  except  the  water  vapor  in  the  air. 

Based  upon  the  above  tests,  it  was  concluded  that  there  would  be 
no  variation  in  the  moisture  content  of  the  sand  during  the  course  of 
the  experiment. 

Direct  Shear  Test 

In  a  cohesionless  soil,  the  angle  of  internal  friction  is  the 
soil  parameter  which  determines  the  shear  resistance.  In  order  to 
determine  0,  a  strain-controlled  direct  shear  test  was  conducted  on  a 
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dense  sample  of  sand  as  outlined  in  Chapter  X  of  Reference  31*  In. the 
direct  shear  test,  a  portion  of  the  sand  sample  is  caused  to  slide  in 
relation  to  the  rest  of  the  saug>le«  The  tests  were  conducted  by  stu¬ 
dents  at  the  University  of  Dayton.  Dayton,  Ohio.  Plotting  the  data 
for  all  the  tests  and  fitting  a  straight  line  to  the  data  results  in 
an  angle  of  internal  friction  of  37° •  This  result  is  shown  in  Fig.  52 
where  each  data  point  is  the  state  of  stress  on  the  soil  sample  at  the 
time  of  shear  failure.  Because  a  constant  area  of  the  sample  is  as¬ 
sumed,  the  scale  is  in  terms  of  pounds  of  force.  To  compute  the  stress 
in  psf,  each  value  must  be  divided  by  1/30  square  foot. 

Mention  should  be  made  of  the  fact  that  the  line,  which  is  really 
a  plot  of  the  Coulomb  shear  formula 

s  *  c  +  p  tan  i  (lh) 

indicates  that  the  sand  apparently  possesses  some  cohesion  as  indi¬ 
cated  by  the  intercept  on  the  axis.  This  is  due  to  the  presence  of 
friction  between  the  tc-p  of  the  shear  box  and  the  bottom  of  the  shear 
box  (Ref  19) . 

Summary 

The  sieve  analysis  indicated  that  the  soil  was  uniformly  graded. 
The  intensification  tests  showed  that  the  variation  in  density  under  the 
conditions  of  the  test  would  be  negligible.  Also,  there  was  negligible 
moisture  in  the  sand  under  the  conditions  of  the  test.  The  direct 
shear  test  indicated  that  the  soil  was  cohesionless  and  had  an  angle 
of  infernal  friction  of  37° •  Based  upon  these  characteristics,  the 
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medium  could  be  considered  to  be  an  ideal  medium  in  which  to  invest 
gate  the  effect  of  gravity  on  crater  formation. 
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Appendix  C 

Tabulated  Data 

This  appendix  contains  a  listing  of  the  dimensions  of  the  craters 
formed  during  the  calibration  and  test  phases  of  the  experiment.  Also 
included  are  the  means  and  variances  of  the  craters  formed  under  the 
various  calibration  and  test  conditions.  All  dimensions  and  depths  of 
burst  listed  in  the  tables  are  expressed  in  inches. 
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TABLE  II 

DIMENSIONS  OF  CRATERS  FORMED 
UNDER  .17  0  CONDITIONS 


0.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

371 

11.750 

2.0C0 

372 

12.000 

2.125 

373 

12.500 

2.250 

37U 

10.250 

1.500 

375 

9.000 

1.250 

376 

9.250 

1.625 

377 

11.250 

2.000 

378 

9-500 

1.375 

379 

10.250 

2.000 

380 

10. 7*0 

1.250 

1.0 

Inch  DOB 


121 

16.250 

2.750 

122 

17.500 

2.750 

123 

17.250 

2.375 

121; 

17.000 

2.125 

125 

18.250 

3.000 

126 

16.250 

2.750 

127 

17.000 

2.625 

128 

17.750 

2.875 

129 

17.500 

3.125 

130 

17.250 

3-000 

0.5 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

251 

9.500 

1.500 

252 

10.500 

1.875 

253 

10.500 

2.000 

25U 

10.500 

2.000 

255 

10.500 

1.875 

256 

10.250 

1.875 

257 

10.750 

2.125 

258 

10.250 

1.750 

259 

11.000 

2.000 

260 

10.750 

2.000 

1.5 

Inch  DOE 


151 

17.000 

2.375 

152 

18.500 

2.500 

153 

18.000 

3.125 

i5lt 

19.000 

3.125 

155 

18.500 

2.875 

156 

17.750 

2.750 

157 

18.250 

2.625 

158 

19. <000 

3.125 

159 

18.250 

3.250 

160 

17.500 

2.750 
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TABLE  II  (CONTINUED)  DIMENSIONS  OF  CRATERS 


2.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

na 

17.500 

2.375 

11*2 

18.000 

2.125 

lU3 

16.750 

2.000 

lkb 

17.500 

2.375 

1 U5 

17.750 

2.750 

Hi6 

18.000 

3.000 

1L7 

18.000 

2.785 

148 

18.500 

2.750 

1U9 

17.500 

2.750 

150 

18 .000 

2.875 

b.O 

Inch  DOB 


181 

12.750 

1.250 

182 

lii.  500 

1.500 

183 

11;.  000 

1.250 

18U 

13.500 

1.125 

185 

Hi. 000 

1.125 

186 

13.125 

0.625 

167 

13.750 

0.750 

188 

12.750 

0.750 

189 

13. 000' 

1.125 

190 

lli. 000 

1.000 

3.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

161 

15.000 

1.875 

162 

15.250 

2.000 

163 

15.250 

1.000 

I6ii 

16.000 

1.87  5 

165 

16.000 

1.875 

166 

16.750 

2.250 

167 

16.250 

1.875 

168 

16.250 

2.125 

169 

15.750 

1.875 

170 

15.750 

2.125 

5.0 

Inch  DOB 


201 

7.000 

.625 

202 

11.000 

.875 

203 

8.500 

.625 

20li 

9.000 

.500 

205 

7.500 

.625 

206 

11.500 

# 

207 

10.000 

.500 

208 

# 

.250 

209 

iff 

.375 

210 

9.000 

.500 

#  Crater  ..ut  measure  able 
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TABLE  III 

DIMENSIONS  OF  CRATERS  FORMED 
UNDER  .38  0  CONDITIONS 


0.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

381 

12.000 

2.375 

382 

10.750 

1.875 

383 

9.000 

1.375 

381* 

9.250 

1.250 

385 

9.750 

1.625 

386 

10.000 

2.250 

387 

10.500 

2.125 

388 

10.250 

1.750 

389 

10.750 

2.000 

390 

11.000 

2.125 

1.0 

Inch  DOB 


1*01 

11*.  250 

2.375 

1*02 

15.750 

2.500 

1*03 

15.750 

2.250 

1*01* 

11*.  750 

2.375 

1*05 

15.000 

2.375 

1*06 

ll*.250 

2.750 

1*07 

15.500 

2.750 

1*08 

16.250 

3.125 

1*09 

16.500 

?  750 

1*10 

16.250 

2.750 

0.5 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

1*11 

15.750 

2.875 

1*12 

11*.250 

2.750 

1*13 

15.000 

2.750 

)|1)| 

H*.25o 

2.375 

1*15 

U*.5oo 

2.625 

1*16 

11*. 750 

2.250 

1*17 

# 

# 

1*18 

H*.  5oo 

2.625 

1*19 

15.000 

2.875 

1*20 

ll*. 000 

2.625 

1.5 

Inch  DOB 


1*31 

15.750 

3.000 

1*32 

16.000 

3.250 

1*33 

16.000 

2.500 

1*31* 

15.750 

3.000 

1*35 

15.500 

2.750 

1*36 

15.000 

2.625 

1*37 

16.000 

2.750 

1*38 

16.000 

2.750 

1*39 

16.500 

2.875 

1*1*0 

15.500 

2.500 
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TABLE  III  (CONTINUED)  DIMENSIONS  OF  CRATERS 


2.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

1*21 

# 

# 

1*22 

15.000 

2.750 

u23 

15.500 

2.250 

1*21* 

15.500 

2,250 

1*25 

15.750 

2.625 

1*26 

16.500 

2.375 

1*27 

16.250 

2.250 

1*20 

16.000 

2.625 

1*29 

15.750 

2.250 

1*30 

15.500 

2.125 

3.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

1*1*1 

12.750 

1.375 

1*1*2 

11*.  500 

1.750 

1*1*3 

11*.  250 

1.625 

1*1*1* 

li*. 5oo 

1.750 

1*1*5 

15.250 

1.875 

1*1*6 

13.500 

1.375 

1*1*7 

13.750 

1.250 

1*1*8 

11*.  750 

2.125 

1*1*9 

13.750 

1.250 

1*50 

H*. 5oo 

1-750 

l*6l 

1*62 

1*63 

1*61* 

1*65 

1*66 

1*6? 

1*68 

1*69 

1*70 


1*.0 


Inch  DOB 

11.750 

0.875 

11.500 

1.000 

11.750 

1.000 

12.250 

0.625 

12.250 

0.875 

11.750 

1.000 

12.000 

0.875 

12.000 

1.000 

12.000 

1.125 

11.750 

1.250 

5.0 

Inch  DOB 


1*81 

7.250 

1*82 

# 

U83 

# 

1*81* 

1*85 

7.000 

1*86 

7.500 

1*87 

8.000 

1*88 

8.000 

1*89 

# 

1*90 

# 

.625 
•  375 
.375 
.625 
•375 
.625 
.625 
.750 
.375 
.375 


#  Crater  not  measureable 


Mech/GSF-65-35 


TABLE  IV 

DIMENSIONS  OF  CRATERS  FORMED 
UNDER  1.0  0  CONDITIONS 


0.0 

Inch  DOB 

Index 

Diameter 

Depth 

Number 

101 

8.750 

1.875 

102 

9.000 

1.875 

103 

8.500 

1.500 

lQli 

9.250 

1.750 

105 

8.250 

1.500 

106 

9.000 

1.500 

107 

8.250 

1.375 

ioe 

8.875 

1.625 

109 

7.375 

1.500 

no 

8.500 

1.500 

in 

9.000 

1.750 

112 

9.000 

1.875 

113 

5.000 

1.750 

llii 

9.250 

1.750 

115 

9.250 

1.625 

116 

9.000 

1.750 

117 

9.500 

1.625 

118 

9.250 

1.875 

119 

9.250 

1.625 

120 

8.750 

1.750 

0.5 


Inch  DOB 

Index 

Diameter 

Depth 

Number 

501 

10.750 

1.875 

502 

11.250 

2.250 

503 

11.750 

2.250 

504 

12.750 

2.250 

505 

11.750 

2.250 

506 

11.000 

2.250 

50? 

12.750 

2.250 

508 

12.000 

2.375 

509 

12.000 

2.250 

510 

12.000 

2.250 

511 

12.500 

2.500 

512 

12.000 

2.125 

513 

12.500 

2.625 

514 

11.750 

2.2^0 

515 

12.250 

2.150 

516 

12.500 

2.125 

517 

12.250 

2.250 

518 

12.500 

2-375 

519 

1Z.000 

2.125 

520 

12.500 

2.250 

Mech/GSF-65-35 


tahle  IV  (continued)  dimensions  of  craters 


Index 

Number 


1.0 

Inch  DOB 


1.5 

Inch  DOB 


Diameter  Depth  index 

Numbar 


Dia*ne  te  r  Depth 


1 

2 

3 

k 

5 

6 

7 

8 
9 

10 

11 

12 

13 

lb 

15 

16 

17 

18 
15 
20 


13.750 

13.375 

12.875 

1U.125 

13.125 

13.125 

13.250 

13-375 

13.125 

13.125 

13.125 

13.500 

13.250 

Hi.  000 

13.500 

U*.000 

13.750 

13.625 

13.50C 

13.750 


2.500  521 
2.625  522 
2.500  523 
2.500  52U 
2.503  525 

2.500  526 
2.375  52V 
2.500  523 
2.375  >29 
2.375  530 
2.250  531 

7.500  53? 
2-250  533 
2.750  534 
2.375  535 

2.500  536 
2.500  537 
2-375  538 
2.750  53} 

2.625  5Uo 


13.25c 

13.500 

13.500 

13.750 
l*i.0Q3 
13.75C 
13. 500 
13.750 
13.750 
13.500 
la. 000 
^4. 000 

1*4.000 

lii.ooa 

13.750 

1*4.250 
1*4.250 
ia.ooo 
Hi. 000 
Hi. 000 


2.500 

2.250 

2.750 

2.375 

2.500 

2.125 

2.250 

2.250 

2.250 

2.500 

2.375 

2.375 

2.250 

2.375 

2.250 

2.250 

2.250 

2.125 

2.250 

2.500 
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TABLE  IV  (CONTINUED)  DBGNSIONS  OF  CRATERS 


2.0 

3-0 

Inch  DOB 

Inch  DOB 

Index 

Diameter 

Depth 

Index 

Diameter 

Depth 

Number 

Number 

21 

13-375 

1.875 

ia 

12.125 

1.625 

22 

13.500 

2.125 

U2 

12.000 

1.125 

23 

13.875 

2.250 

1*3 

11.875 

1.625 

2k 

13.375 

2.000 

Uli 

11.875 

1.125 

25 

13-500 

1.875 

i*S 

12.000 

1.250 

26 

lii.125 

2,000 

U6 

12.000 

1.250 

2? 

13.625 

1.625 

hi 

12.000 

1.375 

26 

13.375 

1.875 

1*8 

11.000 

1.000 

29 

13*125 

1.875 

1*9 

12.25c 

1.375 

30 

13.625 

1.875 

50 

11.625 

1.000 

31 

13-750 

2.125 

51 

13.000 

1.625 

32 

13.375 

1.S75 

52 

11.000 

0.875 

33 

13.875 

2.375 

53 

11.750 

1.375 

3u 

13.000 

1.375 

51* 

12.750 

1.50C 

35 

13.625 

1.875 

55 

12.500 

1.3?5 

36 

13.375 

1.750 

56 

12.625 

2.375 

3? 

13.750 

1.875 

57 

12.750 

1.500 

38 

13.125 

1.750 

58 

11.500 

1.250 

39 

13.750 

2.000 

59 

12.000 

1.125 

uO 

13-37^ 

1.750 

60 

12.125 

1.250 

Kecry  3S?-i5~  35 


TABLE  IV  (CONTINUED)  DIMENSIONS  OF  CRATERS 


U.O 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

61 

10.250 

1.000 

62 

10.250 

0.750 

63 

10.875 

0.750 

6U 

9.8?5 

0.500 

65 

11.000 

1.000 

66 

11.375 

1.000 

67 

10.625 

0.750 

68 

8.500 

0.875 

69 

10.250 

0.875 

70 

10.500 

0.750 

71 

10.000 

1.000 

72 

10.500 

0.750 

73 

10.500 

1.000 

7ii 

10.250 

0.750 

IS 

10.000 

1.000 

76 

10.500 

0.750 

77 

10.000 

0.625 

78 

10.500 

0.750 

79 

9.750 

1.000 

8o 

10.000 

1.000 

11 h 
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TABLE  V 

DIMENSIONS  OF  CRATERS  FORMED 
UNDER  2.5  G  CONDITIONS 


0.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

391 

8.750 

1.500 

392 

8.750 

1.375 

393 

8.000 

1.500 

391+ 

8.500 

1.500 

395 

7.750 

1.625 

396 

8.500 

1.375 

397 

7-750 

1.500 

398 

7.500 

1.375 

399 

7.750 

1.375 

1+00 

7.250 

1.375 

1.0 

Inch  DOB 


21+1 

12.500 

2.375 

2U2 

12.250 

2.000 

21+3 

11.750 

2.000 

21+1+ 

12.000 

2.250 

21+5 

12.250 

2.125 

21+6 

11.750 

2.125 

21+7 

11.750 

2.125 

21+8 

11.250 

2.250 

21+9 

11.750 

2.125 

250 

12.250 

2.375 

0.5 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

251 

9.500 

1.500 

252 

10.500 

1.875 

253 

10.500 

2.000 

251+ 

10.500 

2.000 

255 

10.500 

1.875 

256 

10.250 

1.875 

257 

10.750 

2.125 

258 

10.250 

1.750 

259 

11.000 

2.000 

260 

10.750 

2.000 

1.5 

Inch  DOB 


271 

12.250 

1.875 

272 

12.500 

1.875 

273 

12.250 

1.875 

271+ 

12.250 

1.875 

275 

12.500 

1.750 

276 

12.750 

2.000 

277 

12.750 

2.125 

278 

12.000 

2.000 

279 

13.000 

2.125 

280 

12.500 

1.875 
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TABLE  V  (CONTINUED)  DIMENSIONS  OF  CRATERS 


2.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

26i 

12.000 

1.875 

262 

12.250 

1.750 

263 

12.000 

1.500 

26Ii 

12.500 

1.625 

265 

12.000 

1.625 

266 

12.500 

1.375 

267 

12.000 

i.5oo 

268 

12.500 

1.375 

269 

12.250 

1.500 

270 

12.000 

1.375 

h.o 

Inch  DOB 


301 

8.250 

•750 

302 

8.000 

•  375 

303 

8.750 

.500 

30U 

8.500 

.500 

305 

10.250 

.875 

306 

8.500 

.5oo 

307 

9.250 

.625 

308 

10.000 

.500 

309 

8.500 

.500 

310 

8.750 

.625 

3.0 

Inch  DOB 


Index 

Number 

Diameter 

Depth 

281 

10.250 

1.125 

282 

11.250 

1.000 

283 

10.750 

1.000 

28U 

11.250 

1.125 

285 

10.000 

1.000 

286 

10.750 

1.125 

287 

11.000 

1.000 

288 

10.750 

1.125 

289 

10.000 

1.000 

290 

10.500 

0.875 

K*ch/08T-6S-35 
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TABLE  VII 

DIMENSIONS  OF  CRATERS  FORKED 
DURING  CALIBRATION  TESTS  AT  1  g 

Craters  Formed  In  Thirteen  Inches  of  Sand 


DOB 

Number 

Diameter 

Mean 

Variance 

Depth 

Mean 

Variance 

13-1 

10.625 

0.750 

13-2 

9.500 

0.750 

U.O 

13-  3 

10.250 

10.075 

.160 

1,000 

.350 

.009 

13-U 

10.25C 

0.675 

13-5 

9.750 

0.675 

Craters  Formed  Without  Reference  Wire 


-1-1  13.250  2.500 

-1-2  13.375  2.500 

-1-3  13.625  13.400  .021  2.500  2.550  .010 

-1-li.  13.500  2.50C 

-1-5  13.250  2.750 


1,000 

0  S75 

.635  0*875  -850  .009 

o,75o 

0.750 


k~l 
■h-2 
•4-3 
■U-U 
WW-U-S 


10,000 

9.750 

10.000 

9.500 

11.750 


10.200 
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TABLE  VII  (CONTINUED)  DIMENSIONS  OF  CRATERS  FOR  CALIBRATION 


Craters  Formed  For  Vibration  Teats 


DOB 

Number 

Diameter 

Kean  Variance 

Depth 

Mean 

Variance 

V-l-1 

13.750 

2.125 

V-l-2 

Ik. 000 

2.125 

1.0 

V-l-3 

13.500 

13.U50  .160 

2.125 

2.225 

.015 

V-l-U 

13-000 

2.375 

7-1-5 

13.000 

2.375 

7-4-1 

11.000 

.500 

V-ii-2 

9.5oo 

.625 

fc.O  V-U-3 

9.750 

9.7,50 

.1*50 

.750  .500 

.031 

V-h-U 

9.500 

7-4-5 

9.000 

.250 

11? 


H*ce/UHT-G>-  yz 


TABLE  VII  (CONTINUED)  DBffllSIONS  OF  CRATERS  FOR  CALIBRATION 
Craters  Formed  to  Determine  Effect  of  DOB  Placement 


DOB  Number 

Diameter 

Mean 

Variance  Depth 

Mean 

Variance 

E-l 

10.750 

2.250 

0.5'“'B-3 

11.250 

2.000 

11.750 

11.250 

.100  2.125 

2. 050 

.016 

E-U 

11.250 

1.875 

&-5 

11.250 

2.000 

JS-6  12.000  2.250 

.  .E-7  12.500  2.250 

0.5^J&-8  12.250  12  .LOO  .065  2.250  2.275  .003 

K-9  12.500  2.250 

E-10  12.750  2.375 
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Appendix  D 


This  appendix  contains  photographs  and  profiles  of  typical 
craters  formed  at  each  depth  of  burst  and  in  each  gravity  field.  The 
shape  of  the  profile  was  obtained  from  direct  observation  of  the 
craters,  examination  of  the  cratering  films,  and  the  measured  depth 
and  diameter.  The  crater  rim-to-rim  diameter  and  crater  depth  are 
accurately  represented  while  the  remaining  dimensions  such  as  rim 
width  and  height  are  estimated.  However,  the  photographs  and  profiles 
do  give  an  indication  of  the  tre..d  in  the  crater  size  and  shape  as 
gravity  end  depth  of  burst  are  varied. 

In  certain  of  the  crater  profiles,  a  hummocky  rim  is  noted. 

During  the  experiment,  some  of  the  fallback  would  occasionally  trap 
the  expanding  gases  before  they  vented  into  the  air.  This  resulted 
in  a  secondary  scouring  action  as  the  gases  escaped  from  the  additional 
entrapment.  This  secondary  scouring,  which  appeared  to  be  independent, 
of  gravity,  is  believed  to  be  the  roason  for  the  hummocky  rims. 


121 


Profile  of  Crater  Formed  at  C.C  I.  ch  BCB  at  A 


Profile  of  Crater  Formed  at  0.0  Inch  DOB  and  1*0  g 


Hech/GSF-65-3? 


Figure  57 

Profile  of  Crater  Formed  at  0.5  Inch  DOB  and  0.17  g 


t  0.5  Inch  DOB  and  1.0  g 
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Profile  of  Crater  Formed  at  1,0  Ii-ch  DOB  and  0,33 
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Figure  6.3 

Profile  of  Crater  Formed  at  1.0  inch  DOB  and  1.0 


Profile  of  Crater  Formed  at  1*5  I- eh  DOB  and  0*17  g 


Figure  67 

Profile  cf  Crater  Formed  at  1.5  Inch  DOB  and  1.0  g 


Meeh/GSP-65-3 


f 


Mech/GSF-65-35 


Profile  of  Crater  Formed  at  2.0  l:.ch  DC 3 
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Profile  of  Crater  Formed  at  2.0  Inch  DOB  and  1.0 
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Profile  of  Crater  Forr.oJ  at  3-C  li.ci*  DO!  a..a  C.17  g 
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Figure  76 

Profile  of  Crater  Formed  at  3-0  I -oh  DOB  and 


Profile  of  Crater  Formed  at  3«0  Inch  DOB  and  1.0  g 
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0  Inch  DOB  and  0.3 


Profile  of  Crater  Formed  at  4.0  Inch  DOB  ar*d  0.17  g 
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Profile  of  Crater  Formed  at  4.6  Inch  DOB  and  1.0  g 


Profile  of  Crp*.*»r  Formed  at  5.0  Inch  DOB  and  0.33 
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Figure  83 

Profile  of  Crater  Formed  at  5*0  Inch  DOB  ar.d  1.0  g 


Appendix  E 


Statistical  Analysis  and  Computations 

In  this  section  the  statistical  tests  which  were  used  in  analyz¬ 
ing  the  experimental  data  are  discussed.  The  Kolmogorov- Smirnov  test, 
which  io  disoussed  first,  was  necessary  to  determine  if  the  sample  data 
could  be  assumed  to  be  normally  distributed.  A  normal  distribution 
was  necessary  as  a  requirement  for  using  the  two-sample  t-test  to  de¬ 
termine  if  gravity  had  a  significant  effect  on  crater  dimensions .  The 
results  of  these  two  tests  are  tabulated  in  this  section.  Finally,  an 
explanat;  on  is  given  of  the  least-squares  method  used  in  determining 
the  relationship  between  crater  dimensions  and  gravity. 

Kolmogorov-Smirnov  Test 

Whenever  it  is  desired  to  subject  data  to  some  type  of  statisti¬ 
cal  analysis,  one  of  the  first  problems  is  to  determine  what  proba¬ 
bility  distribution  function  F(x)  the  sample  cumulative  distribution 
function  S(x)  approximates.  If  the  assumption  is  made  that  F(x)  is 
continuous,  a  test  which  may  then  be  used  is  the  Kolmogorov-Smirnov 
(Re:  60:$0). 

In  this  test,  the  maximum  absolute  deviation  Dn  of  the  sample 
S '  )  from  what  is  assumed  to  be  the  true  distribution  function  F(x) 
is  determined  either  by  tabulating  the  two  functions  or  by  graphical 
means.  A  theorem  can  be  proved  which  states  that  as  long  as  F(x)  is 
continuous,  S(x)  will  approximate  it  in  the  same  way.  .low  well  S(x) 
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approximates  F(x)  can  be  determined  by  using  as  a  criterion 
(Ref  60  j  pi) .  As  a  consequence,  Dq  is  a  random  variable  which  can  be 
used  for  constructing  a  confidence  band  for  F(x) .  In  order  to  use 
the  random  variable  Dn,  it  is  necessary  to  know  its  distribution. 
This  can  be  worked  out  for  a  particular  set  of  n  samples  and  desired 
level  of  significance  Od  •  The  tern  level  of  significance  or  a  is 
defined  as  the  probability  of  rejecting  the  stated  statistical  hy¬ 
pothesis  when  it  should  have  been  accepted.  Certain  critical  values 
of  Dq  are  tabulated  in  statistical  texts.  These  critical  values  are 
denoted  by  the  symbol 

t*<;n 

where OC  ■  level  of  significance 

n  -  number  of  observations  In  the  sample 
Therefore,  for  a  given  sample  of  size  n  and  desired  level  of  signifi¬ 
cance,  a  value  of  t^^  can  be  obtained.  For  the  Kolmogorov-Smirnov 
Test  where  is  the  criterion  for  the  fit  of  S(x)  to  F(x),  the  hy¬ 
pothesis  is  made  that 

P  (Dn^t^p  J  F(x)  is  the  true  distribution  function  of  x)  -  OC 

where  P  denotes  the  probability  and  all  the  other  terms  are  as  pre¬ 
viously  defined.  Then  if 

the  hypothesis  is  rejected;  If 
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the  hypothesis  is  accepted  (Ref  60:£l). 

For  the  crater  dimensions  obtained  during  the  experimental  phase, 
the  hypothesis  was  made  that  the  sample  data  approximated  a  normal  or 
Gaussian  distribution.  A  level  of  significance  of  0.01  was  chosen  as 
the  minimum  level  at  which  the  hypothesis  would  be  accepted.  Although 
an  cx  of  0.01  appears  to  be  a  weak  requirement,  the  test  is  so  devised 
that  it  is  easier  to  reject  the  hypothesis  than  to  accept  it  (Ref  60:5l)i 
Therefore,  the  ability  to  accept  the  hypothesis  at  even  a  0.01  level 
of  significance  is  meaningful.  As  will  be  shown  in  the  tabulated  re¬ 
sults  of  the  test,  it  was  possible  to  accept  the  hypothesis  with  an  0( 
of  0.20  in  the  majority  of  the  tests. 

To  illustrate  the  preceding  technique,  a  sample  calculation  is 
done  on  the  diameters  of  craters  formed  under  0.17  g  conditions  at 
1.0  inch  DOB.  These  measurements  are 

16.250,  17-500,  17.250,  17.000,  18.250 

16.250,  17.000,  17.750,  17.500,  17.250 

The  mean  of  this  data  is  17.2  and  the  variance  is  O.3U8.  From  this 
data,  the  cumulative  sample  distribution  can  be  tabulated. 

For  each  diameter,  it  is  necessary  to  compute  F(x)  by  means  of 
the  expression 

F|x) 
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where 


*  .  X-At 

1 "  ~6~  ™ 

and  X-  measured  dimension 
(A=  mean  of  sample  data 
( J  =  standard  deviation  of  sample  data 
The  expression  for  F(x)  has  been  extensively  tabulated.  There¬ 
fore,  by  converting  each  of  the  sample  dimensions  to  a  value  of  t  by 
means  of  Eq.  (16),  the  value  of  F(x)  can  be  obtained  from  tables. 
These  values  along  with  the  cumulative  sample  distribution  S(x)  in 
table  form  are 


Range  of  x 

S(x) 

X 

F(x) 

x-16.250 

0.00 

15.000 

.00 

l6.250ix-17.000 

0.20 

16.250 

.06 

17  • 000:2  x-«  17. 250 

o.Uo 

17.000 

.37 

17-25021x-«17.500 

0.60 

17.250 

.53 

17.500-x-17.750 

0.80 

17.500 

.70 

17.750ix-l8.250 

0.90 

17.750 

.82 

l8.250ix-oo 

1.00 

18.250 

.96 

Since  F(x)  is  an  increasing  continuous  function  and  S(x)  is  an 
increasing  step  function,  the  only  differences  which  must  be  examined 
to  determine  Dn  are  the  differences  that  occur  at  the  discontinuities 
in  S(x).  These  differences  or  deviations  are  seen  in  Fig.  85  which 
shows  a  plot  of  F(x)  and  S(x).  These  deviations,  which  can  be  obtained 
from  the  figure  or  the  table,  are  tabulated  below. 
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Value  of  X 


Left  Hand  End  Point 


Right  Hand  End  Point 


15.000 

.00 

.00 

16.250 

.11* 

-.06 

17.000 

.03 

-.17 

17.250 

.07 

-.13 

17.500 

.10 

-.10 

17.750 

.08 

-.02 

18.250 

.06 

-.01* 

Examination  of  the  differences  gives  the  value  of 

D10  -  Max  |  S(x)  -  F(x)|  -  0.17 

Referring  to  a  table  which  lists  the  values  of  t^  .n  where  n  *  10,  it 
is  found  that 

Pd^O  »  .323)  -  0.20 

This  mans  that  the  hypothesis  of  a  normal  distribution  with  re an  17.2 
and  variance  0.31*8  for  the  sample  data  can  be  accepted  on  a  signifi¬ 
cance  level  of  0.20.  The  ability  to  accept  on  such  a  high  level  of 
significance  is  a  very  positive  indication  that  the  sample  data  is 
normally  distributed. 

The  Kolmogorov-Snirnov  Test  was  performed  on  all  the  crater  di¬ 
mensions.  In  each  case,  the  hypothesis  was  made  that  the  true  dis¬ 
tribution  of  the  sample  data  was  the  normal  distribution.  In  each 
test,  the  highest  possible  (X  wee  used.  However,  in  no  instance  wee 
an  (X  leas  than  0.01  considered  acceptable.  The  results  of  these 
tests  are  tabulated  in  this  section.  From  these  tables  it  can  be  seen 
that  63  out  of  70  hypotheses  were  able  to  be  accepted  with  an  Q!  of 


0.20.  In  only  one  instance  did  an  OC  of  O.GJ,  nave  to  be  as^d .  In 
no  instance  did  the  hypothesis  have  to  be  rejected.  Therefore,  the 
data  is  normally  distributed. 
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TABLE  VI  IT 
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RESULTS  OF  KOLMOOOROV-SMIRNOV 
TESTS  ON  CRATER  DIMENSIONS 

Diameter 


Value 

DOB 

Deviation 

Critical 

Level  c 

of 

in 

of  Sample 

Value  of 

Signific; 

Gravity 

Inches 

Bn 

a 

0.0 

•111 

.323 

,20 

0.5 

.25 

•323 

.20 

1.0 

.17 

.323 

.20 

-1? 

1.5 

.15 

.323 

.20 

2.0 

.21 

.323 

.20 

3.0 

.19 

.323 

.20 

L.o 

.19 

.323 

.20 

5.0 

.18 

*359 

.20 

0.0 

.11 

.323 

.20 

0.5 

.19 

.323 

.20 

1.0 

.17 

.323 

.20 

,36 

1.5 

.20 

.323 

.20 

2.0 

.17 

.339 

.20 

3.0 

.20 

.323 

.20 

li.O 

.2k 

.323 

.2G 

5.0 

.27 

M7 

.20 

O.Q 

.22 

.232 

.20 

0.5 

.17 

.232 

.20 

1.0 

•lh 

.232 

.20 

1.0 

1.5 

.25 

.232 

.20 

2.0 

.16 

.232 

.20 

3.0 

.13 

.232 

.20 

H.o 

.17 

.232 

.20 

0.0 

.22 

.323 

.20 

o.5 

.25 

.323 

.20 

1.0 

.22 

.323 

.20 

2.5 

1.5 

.19 

.323 

.20 

2.0 

•  32 

.323 

.20 

3.0 

.19 

.323 

.20 

h.O 

.27 

.323 

.20 
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TABLE  Till  (CONTINUED)  RESULTS  OF  TESTS 


Depth 

Value 

DOB 

Deviation 

Critical 

Level  of 

of 

in 

of  Sample 

Value  of 

Significance 

Gravity 

Inches 

Dn  > 

a 

0.0 

.2? 

.323 

.20 

0.5 

.22 

.323 

.20 

1.0 

.21 

.323 

.20 

.17 

1.5 

.26 

.323 

.20 

2.0 

.29 

.323 

.20 

3.0 

•  39 

.U09 

.05 

ii.O 

.22 

.323 

.20 

5.0 

.20 

.339 

.20 

0.0 

.16 

.323 

.20 

0.5 

.15 

.323 

.20 

1.0 

.21 

.323 

.20 

.38 

1.5 

.19 

.323 

.20 

2.0 

.17 

.339 

.20 

3.0 

.21 

.323 

.20 

U.o 

.21 

.323 

.20 

5.0 

,6U 

.669 

.01 

0.0 

.20 

.232 

.20 

o.5 

.32 

.329 

.02 

1.0 

.25 

.265 

.10 

1.0 

1.5 

.28 

•  29U 

.05 

2.0 

.21 

.232 

.20 

3.0 

.17 

.232 

.20 

U.o 

.26 

.265 

.10 

0.0 

.32 

.323 

.20 

o.5 

•  2U 

.323 

.20 

X  .0 

.25 

.323 

.20 

2.5 

1.5 

.19 

•  323 

.20 

2.0 

.2U 

.323 

.20 

3.0 

.28 

.323 

.20 

U.o 

.31 

.323 

.20 
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TABLE  IX 

RESULTS  OF  KOLMOGORGV-SMI RNO V  TESTS  ON 
CRATERS  USED  FOR  CALIBRATION  PURPOSES 


Condition 

DOB 

Dimension 

13  Inches 

U.o 

Diameter 

of  Sand 

U.o 

Depth 

Without 

1.0 

Diameter 

Reference 

1.0 

Depth 

Wire 

U.o 

Diameter 

u.o 

Depth 

Aircraft- 

1.0 

Diameter 

Titration 

1.0 

Depth 

U.o 

Diameter 

U.o 

Depth 

W** 

ex# 

.27 

.UU7 

.20 

,26 

•UU7 

.20 

,25 

.UU7 

.20 

,U9 

.519 

.05 

.Uo 

.UU7 

.20 

.26 

•UU7 

.20 

.27 

•UU7 

.20 

.39 

.UU7 

.20 

•  30 

•UU7 

.20 

.16 

•UU7 

.20 

•w-  Deviation  of  sample 
##  Critical  value  of  deviation 

#  Level  of  significance 
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t-Test 


In  order  to  justifiably  conclude  that  gravity  did  or  did  not 
affect  the  crater  dimensions,  it  was  decided  to  resort  to  a  statisti¬ 
cal  test.  The  desired  test  would  be  one  which  would  compare  the  mean 
diameter  or  depth  of  craters  formed  at  the  same  DOB  but  under  differ¬ 
ent  gravity  conditions.  If  the  test  showed  that  there  was  a  signifi¬ 
cant  difference  between  the  two  means,  then  it  could  be  stated  with  a 
certain  probability  or  level  of  significance  that  gravity  had  an  effect 
on  the  crater  dimension  at  that  DOB.  Conversely,  if  the  test  resulted 
in  the  probability  that  the  means  were  equal,  then  gravity  did  not 
affect  the  crater  dimensions. 

Knowing  that  the  sample  data  was  normally  distributed,  it  was 
possible  to  use  the  two-sample  t-test.  An  advantage  of  the  t-test 
is  that  the  sample  populations  may  be  small.  However,  in  order  to  use 
the  two- sample  t-test,  it  is  necessary  to  make  the  assumption  that  the 
two  distributions  being  compared  have  equal  standard  deviations,  i.e., 
(jf  —  (j  2  •  Looking  at  the  sample  variances,  it  is  seen  that  in  only 
one  instance  do  the  variances  differ  by  more  than  a  factor  of  ten. 
Additionally,  if  the  sample  sizes  of  the  two  populations  are  not 
markedly  different,  then  the  test  is  not  as  sensitive  to  the  assumption 
that  the  sample  variances  are  equal  (Ref  8 j  239) .  Therefore,  the  use 
of  the  t-test  is  appropriate. 

The  t-test  is  derived  by  means  of  the  likelihood  ratio  technique 
which  yields  a  ratio  based  on  the  t  distribution  with  n^  +  n£  -  2 
degrees  of  freedom  and  the  statistic 
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fVI|  sf » (TVyl)  4  OT 

n|+n2-2  V"i  n2 

where  x  =  means  of  sample  data 

(3  -  value  by  which  true  means  differ 
n  =  sample  size 

p 

s  =  variance  of  sample  data 

1,2  =  subscripts  which  refer  to  sample  populations  being  compared 
In  this  test,  the  common  procedure  is  to  state  the  null  hypothesis 


H  :  U.-U  =  6 

o  I  2 

against  one  of  the  alternatives 

5,  ^-6;  W5 

A  value  of  t  is  then  calculated  from  Eq.  (17).  This  value  is 
compared  with  the  value  of  the  random  variable  which  has  the  t  dis¬ 
tribution  with  n^  +  np  -  2  degrees  of  freedom  and  at  a  chosen  level 
of  significance  OC  .  Depending  on  which  alternative  was  chosen,  the 
magnitude  of  t^  when  compared  with  the  calculated  t  will  determine 
whether  the  hypothesis  is  accepted  or  rejected. 

For  this  paper,  the  hypothesis  is 

H.^-VO 

where  (li  m  mean  dimension  at  a  given  1X)B  for  the  lower  value  of 
'  gravity 

mean  dimension  at  a  given  DOB  for  the  higher  value  of 
gravity 

against  the  alternative 
The  maximum  acceptable  level  of  significance  was  0.10.  In  the  test,  if 
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then  the  null  hypothesis  must  be  rejected.  If  the  hypothesis  is  re¬ 
jected,  this  signifies  that  gravity  did  effect  the  crater  dimensions 
at  that  LOB.  Remembering  that  QC  is  the  probability  of  rejecting  the 
hypothesis  when  in  fact  it  should  have  been  accepted,  an  OC  =  0.10 
indicates  that  a  0.10  probability  is  to  be  allowed  in  making  the  im¬ 
proper  conclusion  that  gravity  had  an  effect  on  the  dimensions.  How¬ 
ever,  in  every  t-test,  the  lowest OC  possible  was  chosen  since  the 
lower  the  value  of  OC  ,  the  lower  was  the  probability  of  making  the 
incorrect  conclusion  that  gravity  had  an  effect  on  the  crater  dimen¬ 
sions.  To  illustrate  the  technique,  a  sample  calculation  is  performed. 

It  is  desired  to  determine  if  gravity  had  an  effect  on  the 
diameter  of  craters  formed  under  the  following  conditions: 


g 

DOB 

Mean 

Diameter 

Variance 

Sample  Size 

0.17 

ll.O 

13.550 

0-066 

10 

O.38 

U.o 

11.900 

0.025 

10 

Based  on  the  tabled  data,  t  is  calculated  by  means  of  Eq.  (17)  and 
is  8.^09.  Consulting  a  t-distribution  table  and  noting  the  value  of 
t  for  18  degrees  of  freedom,  one  gets 

‘.001  ■  3-61 

Therefore,  since  t  »\ooi*  the  null  hypothesis  must  be  rejected.  The 
conclusion  is  that  there  is  a  difference  in  the  means  and  that  this 
difference  is  due  to  the  different  gravity  fields.  Since  OC  *  0.001, 
there  is  a  0.001  probability  that  this  conclusion  is  incorrect. 
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The  results  of  all  the  two-sample  t-tests  are  listed  in  Table  X. 
All  possible  combinations  between  craters  formed  at  the  same  DOB  but 
in  different  gravity  fields  are  considered. 

The  two-sample  t-test  was  also  used  in  comparing  the  means  of 
the  craters  formed  during  the  calibration  phase  of  the  experiment  (see 
Appendix  A) .  For  all  the  t-tests  which  were  done  on  the  calibration 
data,  the  hypothesis  was 

H:^-VO 

where /J.f=  mean  of  one  group  of  sample  data 

ur  mean  of  second  group  of  sample  data 
The  alternative  was  stated  as 

For  the  calibration  comparisons*  it  was  desired  that  the  stated  hy¬ 
pothesis  be  accepted  whenever  possible.  For  this  reason*  the  minimum 
oc  was  chosen  to  be  0.002.  The  necessary  condition  for  accepting 
the  hypothesis  was 


*  to</2 

The  results  of  the  t-tests  on  the  craters  formed  during  the  calibration 
phase  show  that  the  hypothesis  that  the  means  are  equal  can  be  ac¬ 
cepted  in  all  cases  except  two.  These  two  cases  are  in  the  comparison 
of  the  depths  of  craters  subjected  to  aircraft  vibration.  The  results 
of  the  t-tests  are  tabulated  in  Table  XI . 

Least- Squares  Calculations 

For  the  reasons  already  stated  in  Section  IV*  the  equation  which 
relates  crater  dimensions  at  a  given  DOB  to  the  gravity  field  was 
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TABLE  X 

RESULTS  OF  T-TEST  ON  EFFECT  OF 
GRAVITY  ON  CRATER  DIMENSIONS 


Diameter 


DOB 


0.0 


0.5 


1.0 


1.5 


2.0 


Gravity 

Gravity 

Field 

Field 

.17 

•  38 

.17 

1.00 

.17 

2.50 

.38 

1.00 

.38 

2.50 

1.00 

2.50 

.17 

.38 

•17 

1.00 

.17 

2.50 

•  38 

1.00 

•  38 

2.50 

1.00 

2.50 

•  17 

.38 

.17 

1.00 

.17 

2.50 

.38 

1.00 

•  38 

2.50 

1.00 

2.50 

.17 

CO 

rr\ 

• 

.17 

1.00 

.17 

2.50 

.38 

1.00 

.38 

2.50 

1.00 

2.50 

.17 

.38 

.17 

1.00 

.17 

2.50 

.38 

1.00 

.38 

2.50 

1.00 

2.50 

Sample 

Tabulated 

t 

0.724 

1.33 

6.124 

3.41 

6.551 

3.61 

6.189 

3.41 

7-339 

3.61 

4.244 

3.41 

-1.382 

1.33 

4.587 

3.41 

21.138 

3.61 

12.435 

3.42 

16.908 

3.65 

8.306 

3.41 

4.511 

3.61 

22.103 

3.41 

24.191 

3.61 

8.220 

3.41 

11.586 

3.61 

11.360 

3.41 

10.506 

3.61 

27.777 

3.41 

27.051 

3.61 

16.598 

3.41 

21.972 

3.61 

12.914 

3.41 

9.958 

3.65 

32.049 

3.41 

35.243 

3.61 

14.937 

3-65 

23.234 

3.65 

13.217 

3.41 

Level  of 
Significance 

a 

* 

.001 

.001 

.001 

.001 

.001 

* 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 

.001 
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TABLE  X  (CONTINUED)  RESULTS  OF  T-TEST  ON  GRAVITY 

Diameter 


Level  of 


Gravity- 

Gravity 

Sample 

Tabulated 

Significance 

DOB 

Field 

Field 

t 

** 

(X 

.17 

.38 

5.989 

3.61 

.001 

.17 

1.00 

18.523 

3.U1 

.001 

3.0 

.17 

2.50 

23 .1*09 

3.61 

.001 

.38 

1.00 

9.U93 

3.1a 

.001 

.38 

2.50 

15-909 

3.61 

.001 

1.00 

2.50 

7.292 

3.U1 

.001 

.17 

.38 

8.509 

3.61 

.001 

.17 

1.00 

III-. 957 

3.1a 

.001 

li.O 

.17 

2.50 

16.386 

3.61 

.001 

.38 

1.00 

8.696 

3.1a 

.001 

.38 

2.50 

10.603 

3.61 

.001 

1.00 

2.50 

5-91ii 

3.1a 

.001 

5.0 

.17 

.38 

2.388 

2.20 

.025 

Depth 

.17 

.38 

-0.850 

1.33 

* 

.17 

1.00 

0.7U8 

1.33 

* 

0.0 

.17 

2.50 

3.688 

3.61 

.001 

.38 

1.00 

2.266 

2.05 

.025 

.38 

2.$0 

3.699 

3.61 

.001 

1.00 

2.50 

U.311* 

3.1a 

.001 

.17 

.38 

-5.18 : 

3.65 

* 

.17 

1.00 

-2.965 

2.77 

* 

o.5 

.17 

2.50 

8.962 

3.61 

.001 

.38 

1.00 

5.81*7 

3.U2 

.001 

.38 

2.50 

7.979 

3.65 

.001 

1.00 

2.50 

6.028 

3.U1 

.001 

.17 

.38 

0.651 

1.33 

* 

.17 

1.00 

2.291 

2.06 

.025 

1.0 

.17 

2.50 

3.61*2 

3.61 

.001 

•38 

1.00 

1.355 

1.31 

.10 

.38 

2.50 

2.510 

2.10 

.025 

1.00 

2.50 

5.997 

3.U1 

.001 
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TABLE  I  (CONTINUED)  RESULTS  OF  T-TEST  ON  GRAVITY 

Depth 


Level  of 


Gravity 

Gravity 

Sample 

Tabulated 

Significance 

DOB 

Field 

Field 

t 

V 

oc 

.17 

•  38 

O.U35 

1.33 

* 

.17 

1.00 

6.532 

3-la 

.001 

1.5 

.17 

2.50 

9.1*07 

3.61 

.001 

.38 

1.00 

6.7UO 

3-1*1 

.001 

.38 

2.50 

10.775 

3.61 

.001 

1.00 

2.50 

7.1*72 

3.1*1 

.001 

.17 

.38 

1.567 

1.33 

.10 

.17 

1.00 

6.928 

3-la 

.001 

2.0 

.17 

2.50 

3.881 

3.61 

.001 

.38 

1.00 

8.897 

3-65 

.001 

to 

rr\ 

• 

2.50 

9.636 

3.65 

.001 

1.00 

2. .50 

1*.1*93 

3.1*1 

.001 

.17 

.38 

2.01*1* 

1.73 

.05 

.17 

1.00 

6.313 

3-la 

.001 

3.0 

.17 

2.50 

8.068 

3.61 

.001 

.38 

1.00 

3.720 

3.1a 

.001 

•38 

2.50 

5.1*58 

3.61 

.001 

1.00 

2.r  ' 

3.831 

3.U1 

.001 

.17 

.38 

0.912 

1.33 

* 

.17 

1.00 

2.798 

2.79 

.005 

U.O 

.17 

2.50 

5.152 

3.61 

.001 

•38 

1.00 

2.028 

1.70 

.05 

.38 

2.50 

U.208 

3.61 

.0C1 

1.00 

2.50 

U.789 

3.U1 

.001 

5.0 

.17 

.38 

0.1*09 

1.33 

* 

*  Hypothesis  that  naans  are  equal  vaa  accepted 
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TABLE  H 

RESULTS  OF  T-TEST  ON  CRATERS  FORMED 
DURING  CALIBRATION  PHASE  OF  EXPERIMENT 

Depth  of  Send 


DOB 

Dimension 

Sample 

Tabulated 

Level  of 

t 

V 

Significance 

a 

Diameter 

.7U2 

.86 

.U0 

U.O 

Depth 

.085 

.26 

.80 

Effect  of  Reference  Wire 

1.0 

Diameter 

.399 

.53 

.60 

1.0 

Depth 

1.071 

1.32 

.20 

U.O 

Diameter 

,2U6 

.26 

.80 

U.O 

Depth 

.085 

.26 

.00 

Ef 

feet  of  Aircraft  Vibration 

1.0 

Diameter 

.035 

.26 

.80 

1.0 

Depth 

3.873 

♦ 

* 

U.o 

Diameter 

1.798 

2.07 

.05 

U.o 

Depth 

U.U82 

* 

# 

*  Hypo  thesis  that  wans  are  the  sane  cjirsnot  be  accepted 
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assumed.  to  be 


D  =  Ae-g2  (18) 

where  D  =  dimension  of  the  crater 

A  =  dimension  at  aero  gravity 
g  =  acceleration  due  to  gravity 
B  =  measure  of  gravity  effectiveness 

The  method  of  least-squares  was  used  to  find  the  values  of  A  and 
3  which  would  best  fit  the  experimental  data.  This  meant  that  it  was 
necessary  to  find  the  minimum  value  of 


4 

*  r  B  -j  2 

D  =  )  [  D±  -  Ae"gi  J  (19) 

t-/ 

where  ET*  =  least-squares  deviation 

Taking  the  partial  derivative  cf  Eq.  (19)  with  respect  to  A  and  3  and 
setting  the  result  equal  to  zero  gives  the  expressions 

tM>*-  -  £  2  [  3i  '  Ae'8?  ]  [  e"Si  ]  -  0  '(20} 

^  (.=/  "  . 


dn  = 


4 

I 

L=( 


Di  - 


3  b 

Ae"gi  j  [  e“si  ]  [  g2  ]  [  In  ]  =0  (21) 


Instead  of  solving  these  equations  simultaneously,  it  was  decided  to 
solve  for  A  in  term3  of  B  from  Eq .  ( 20) .  Then,  by  assuming  a  value  for 
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